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ABSTRACT 

We present new infrared observations of the central regions of the starburst galaxy M 82. The 
observations consist of near-infrared integral field spectroscopy in the H- and if-band obtained 
with the MPE 3D instrument, and of A = 2.4 — 45 /xm spectroscopy from the Short Wavelength 
Spectrometer (SWS) on board the Infrared Space Observatory . These measurements are used, 
together with data from the literature, to (1) re-examine the controversial issue of extinction, 
(2) determine the physical conditions of the interstellar medium (ISM) within the star-forming 
regions, and (3) characterize the composition of the stellar populations. Our results provide a 
set of constraints for detailed starburst modeling which we present in a companion paper. 

We find that purely foreground extinction cannot reproduce the global relative intensities of 
H recombination lines from optical to radio wavelengths. A good fit is provided by a homoge- 
neous mixture of dust and sources, and with a visual extinction of Ay = 52 mag. The SWS 
data provide evidence for deviations from commonly assumed extinction laws between 3 firm and 
10 [im. The fine-structure lines of Ne, Ar, and S detected with SWS imply an electron density 
of w 300 cm~ 3 , and abundance ratios Ne/H and Ar/H nearly solar and S/H about one- fourth 
solar. The excitation of the ionized gas indicates an average effective temperature for the OB 
stars of 37400 K, with little spatial variation across the starburst regions. We find that a random 
distribution of closely packed gas clouds and ionizing clusters, and an ionization parameter of 
f» 10~ 2 ' 3 represent well the star-forming regions on spatial scales ranging from a few tens to a 
few hundreds of parsecs. From detailed population synthesis and the mass-to- if -light ratio, we 
conclude that the near-infrared continuum emission across the starburst regions is dominated by 
red supergiants with average effective temperatures ranging from 3600 K to 4500 K, and roughly 
solar metallicity. Our data rule out significant contributions from older, metal-rich giants in the 
central few tens of parsecs of M82. 

Subject headings: galaxies: individual (M 82) — galaxies: ISM — galaxies: starburst — galaxies: stellar 
content — infrared: galaxies 
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1. INTRODUCTION 

M82 is one of the nearest (3.3 Mpc; Freed- 
man & Madore 1988; 1" = 16 pc), brightest, 
and best-studied starburst galaxies. It has long 
been considered as the archetype of this class of 
objects, and has often been used as test labora- 
tory for starburst theories (see e.g. Telesco 1988 
and Rieke et al. 1993 for reviews). The central 
500 pc of M 82 harbours the most active starburst 
regions and has thus received particular attention 
in the past. Most of M82's infrared luminosity of 
3 x 10 10 L originates from this "starburst core," 
which is severely obscured at optical and ultravi- 
olet wavelengths. 

The qualitative picture of M 82 features the fol- 
lowing. A prominent nucleus, a stellar disk, and 
a kiloparsec-long stellar bar are revealed by near- 
infrared observations {e.g. Telesco et al. 1991; 
McLeod et al. 1993; Larkin et al. 1994). The 
molecular gas resides mainly in a rotating ring or 
tightly wound spiral arms and in an inner spiral 
arm at radii of ~ 400 pc and 125 pc, respec- 
tively {e.g. Shen & Lo 1995; Seaquist, Frayer, & 
Bell 1998; Neininger et al. 1998). The H II re- 
gions are concentrated in a smaller rotating ring- 
like structure of radius « 85 pc, and along the 
stellar bar at larger radii (e. g. Larkin et al. 1994; 
Achtermann & Lacy 1995). HST observations 
have resolved over a hundred compact and lumi- 
nous "super star clusters" across the central kilo- 
parsec (O'Connell et al. 1995; de Grijs, O'Connell, 
& Gallagher 2000). An important series of young 
compact radio supernova remnants extends along 
the galactic plane over 600 pc {e.g. Kronberg, 
Biermann, & Schwab 1985; Muxlow et al. 1994; 
Pedlar et al. 1999) and a bipolar outflow traces 
a starburst wind out to several kiloparsecs {e.g. 
Bregman, Schulman, & Tomisaka 1995; Shopbell 
& Bland-Hawthorn 1998; Lehnert, Heckman, & 
Weaver 1999; Cappi et al. 1999). 

The triggering of starburst activity in M 82 is 
generally attributed to tidal interaction with its 

1 Current address: Leiden Observatory, PO Box 9513, 
2300 RA Leiden, The Netherlands 

1 Based on observations with ISO, an ESA project with 
instruments funded by ESA Member States (especially the 
PI countries: France, Germany, the Netherlands, and the 
United Kingdom) and with the participation of ISAS and 
NASA. The SWS is a joint project of SRON and MPE. 



massive neighbour M81 ~ 10 years ago or, alter- 
natively, to the stellar bar which may itself have 
been induced by the interaction {e.g. Gottesman 
& Weliachcw 1977; Lo et al. 1987; Yun, Ho, & 
Lo 1993, 1994; Telesco et al. 1991; Achtermann 
& Lacy 1995). Beginning with the seminal pa- 
per by Rieke et al. (1980), several authors have 
applied evolutionary synthesis modeling to under- 
stand the nature and evolution of starburst activ- 
ity in M82 {e.g. Bernlohr 1992; Rieke et al. 1993; 
Doane & Mathews 1993; Satyapal et al. 1997). 
However, despite extensive studies, crucial issues 
remain open concerning in particular the compo- 
sition of the stellar population and spatial varia- 
tions thereof, the initial mass function of the stars 
formed in the starburst, and the spatial and tem- 
poral evolution of starburst activity. 

In order to address the above issues quantita- 
tively, we have obtained new observations of the 
central regions of M 82 consisting of near-infrared 
H- and if-band integral field spectroscopy with 
the Max-Planck-Institut fur extraterrestrischc 
Physik (MPE) 3D instrument (Weitzel et al. 
1996), and mid- infrared spectroscopy from the 
Short Wavelength Spectrometer (SWS; de Graauw 
et al. 1996) on board the Infrared Space Observa- 
tory {ISO; Kesslcr et al. 1996). The 3D data pro- 
vide detailed information on small spatial scales 
from key features tracing the stellar populations 
and the interstellar medium (ISM) while the SWS 
data cover the entire A = 2.4 — 45 /im range con- 
taining a wealth of additional starburst signatures. 
The 3D and SWS data allow us to apply a variety 
of new and essential diagnostics which we use, to- 
gether with results from the literature, to address 
several controversial issues and reveal additional 
aspects of M82. 

In this paper, we focus on characterizing the 
physical conditions of the ISM (gas density, abun- 
dances, extinction), the composition of the stellar 
populations (hot massive stars and cool evolved 
stars), and the relative distribution of the gas 
clouds and stellar clusters. In a companion pa- 
per (Forster Schreiber et al. 2000, in preparation; 
hereafter paper 2), we will apply recent starburst 
models to our results to constrain quantitatively 
the star formation parameters and the detailed 
starburst history in the central regions of M82. 

The paper is organized as follows. Section 2 de- 
scribes the observations and data reduction pro- 



2 



ccdure and section 3 presents the results. The 
physical conditions of the ISM are derived in sec- 
tion 4 together with the composition of the mas- 
sive star population, while the composition of the 
cool evolved stellar population is investigated in 
section 5. Section 6 summarizes and discusses the 
results of our analysis. 

2. OBSERVATIONS AND DATA RE- 
DUCTION 

2.1. Near-infrared Data 

We observed M 82 at near-infrared wavelengths 
using the MPE integral field spectrometer 3D 
(Weitzel et al. 1996). A first data set was obtained 
at the 3.5 m telescope in Calar Alto, Spain, on 
1995 January 13, 14, 16, and 21. The observations 
were completed at the 4.2 m William-Herschel- 
Telescope in La Palma, Canary Islands, on 1996 
January 6. 3D slices the focal plane into 16 paral- 
lel "slits" which are imaged and dispersed in wave- 
length by a grism onto a 256x256 HgCdTc NIC- 
MOS 3 array, providing simultaneously the en- 
tire H- or if-band spectrum of each spatial pixel. 
Nyquist-sampled spectra are achieved by dither- 
ing the spectral sampling by half a pixel on alter- 
nate data sets, which are interleaved in wavelength 
after the data reduction. For the M82 observa- 
tions, the instrument setup provided a pixel scale 
of 0.5" pixel" 1 with a field of view of 8" x 8", and 
a spectral resolution of R = A/ A A w 1015 and 
R « 830 in the H- and if-band, respectively. 

In order to sample representative starburst re- 
gions in M82, we selected for the observations an 
area approximately 260 pc x 160 pc (correspond- 
ing to 16" x 10"), nearly parallel to the plane of 
the galaxy. The regions mapped include the nu- 
cleus and extend to the west up to the inner edge 
of the molecular ring (see figure 1). The 3D raster 
consists of four different fields, with adjacent fields 
overlapping by « 3". Each field was observed in an 
object- sky- sky- object sequence, with the off-source 
frames taken on blank portions of the sky 2' — 4' 
away in the north and south directions. The to- 
tal on-source integration times per field and wave- 
length channel were 8—15 and 10 minutes for 
the H- and if-band data, respectively. Typical 
single-frame exposures were 60 — 100 s. For at- 
mospheric calibration, we observed B-, late-F, or 
early-G dwarf stars each night, before and/or af- 



ter M82. The seeing varied from 1" to 1.5" during 
the observations. Table 1 summarizes the log of 
the observations. 

We carried out the data reduction using the 
3D data analysis package, within the GIPSY 
environment (van der Hulst et al. 1992). We 
first performed correction for the non-linear sig- 
nal response of the detector, dark current and 
background subtraction, spatial and spectral flat- 
fielding, wavelength calibration, rearrangement of 
the data in a three-dimensional cube, and bad 
pixel correction following standard procedures as 
described by Weitzel et al. (1996). We reduced the 
reference stars data in the same manner. We then 
corrected the M82 data cubes for atmospheric 
transmission by division with calibration spectra 
obtained as follows. 

In the if-band, we first ratioed the stars spec- 
tra with a black-body curve of the appropriate 
temperature given their spectral type. We re- 
moved the intrinsic stellar features by division 
with the normalized spectrum of the G3 V star 
from the Kleinmann & Hall (1986) atlas, con- 
volved to the resolution of 3D. This template has 
similar absorption line strengths as our F5 V and 
GO V calibrators except for the Br7 line, which 
we removed by linear interpolation. In the H- 
band, no appropriate template spectra were digi- 
tally available, so we constructed composite cal- 
ibration spectra as follows. We generated syn- 
thetic transmission spectra between 1.55 /xm and 
1.75 /xm at the proper zenithal distance using the 
program ATRAN (Lord 1992). At the edges of 
the TJ-band, the transmission is very sensitive to 
the actual atmospheric conditions. Since our ref- 
erence stars have no important features bluewards 
of 1.55 /im and redwards of 1.75 /im (within the 
range observed with 3D), we ratioed these portions 
of their spectrum with the appropriate blackbody 
curve and connected them to the synthetic spec- 
tra. 

The residuals from intrinsic features of the ref- 
erence stars are less than 1%. Those from telluric 
features do not exceed 3% — 4% for most of the 
H- and if-band but amount up to 20% between 
1.9 /zm and 2.0 ^m and at both edges of the H- 
band due to spatial and temporal variability of 
the atmospheric transmission. By fine-tuning the 
airmass for the calibration spectra with the help of 
ATRAN, we reduced the spurious features to less 
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than 1% in the if-band and middle of the iJ-band, 
and less than 10% at the edges of the if-band. 

For proper mosaicking, we smoothed the result- 
ing data cubes with a two-dimensional gaussian 
profile to a common spatial resolution of 1.5". The 
small scale structure is thus smeared out in the 
higher resolution fields but no important spatial 
features are lost. We combined the data cubes 
to produce the final mosaic, adjusting the mean 
broad-band flux in the overlapping areas to a com- 
mon level. We performed absolute flux calibra- 
tion based on the broad-band photometric mea- 
surements of Riekc et al. (1980). We estimate the 
uncertainties on the absolute fluxes to be 15% in 
the H-band and 10% in the if-band. They include 
possible systematic errors in the relative spatial 
and spectral flux distribution which may occur in 
the background subtraction, in the correction for 
telluric absorption, and in the mosaicking (see also 
section 3.2). 

2.2. Mid-infrared Data 

We obtained the mid-infrared spectrum of M 82 
with the SWS (de Graauw et al. 1996) on board 
ISO (Kessler et al. 1996) during revolution 116 
on 1996 March 12. The grating scan mode AOT 
SWS01 was used to cover the entire SWS range 
from 2.4 fim to 45 fim. The slowest scan speed 
was selected to obtain the highest spectral resolu- 
tion possible in this mode (R <~ 500 — 1000). In 
addition, individual lines were scanned with the 
grating line profile mode AOT SWS02 to get full 
resolving power (R <~ 1000 — 2000) and improved 
sensitivity for the key lines used in the data inter- 
pretation. 

The SWS rectangular aperture was centered 
on the western mid-infrared emission peak (ac- 
tual pointing at ai 950 : 09 h 51 m 42!20, 5 1950 - 
+69°55'00'.'74). The major axis of the aperture 
was oriented at a position angle of 64? 5, nearly 
parallel to the galactic plane of M 82. The aper- 
ture size varies from 14" x 20" at short wavelengths 
to 20" x 33" at long wavelengths. The SWS field of 
view thus includes the regions mapped with 3D, 
extending out to a maximum distance of about 
350 pc from the nucleus (see figure 1). The SWS 
full scan took in total l h 50 m . The on-source in- 
tegration time for the individual line scans varied 
between 100 s for most lines and 600 s for the 
weakest lines. 



We reduced the data with the SWS Interac- 
tive Analysis package (SIA), using special inter- 
active extensions for glitch tail removal, dark sub- 
traction, up-down correction, flat-fielding, and de- 
fringing. As the combined spectra of the 12 de- 
tectors of each SWS band are oversampled, we re- 
binned the final data to the proper instrumental 
resolution. We performed the calibration using the 
calibration tables as of 1998 February 15, equiva- 
lent to OLP version 7.0. The overall accuracy of 
the line and continuum fluxes is estimated to be 
w 10% - 20% (Schacidt et al. 1996; Salama et al. 
1997). 

3. RESULTS 

3.1. Near-infrared Spectra 

We selected three individual regions within the 
3D field of view for a detailed analysis: the nu- 
cleus, as well as two regions covering the Br7 
sources located approximately 10" and 5" south- 
west from the nucleus and which we designate as 
M82:Brl and M82:Br2, respectively For conve- 
nience, we will refer to the latter regions as "Bl" 
and "B2" in the rest of this paper. Bl is close 
to the molecular ring and B2 coincides with the 
western mid-infrared emission peak. These three 
positions were chosen because they sample, along 
the galactic plane of M82, representative regions 
with different spectral properties indicative of dif- 
ferent stellar populations as described below. 

We extracted the H- and JT-band spectra of the 
individual regions from the 3D data cubes using 
synthetic apertures of 2.25" x 2.25", correspond- 
ing to 35 pc x 35 pc. They are shown in figure 2 
together with those for the entire 3D field of view. 
All spectra exhibit the typical signatures of star- 
burst activity, including H and He recombination 
lines tracing photoionized nebulae, ro-vibrational 
lines from warm H 2 , [Fe II] transition lines orig- 
inating primarily in iron-enriched shocked mate- 
rial, and numerous atomic and molecular absorp- 
tion features produced in the atmosphere of cool 
stars. However, the intensity of the features rel- 
ative to the continuum varies with position. The 
emission lines become stronger along the sequence 
Nucleus — > B2 — > Bl while the stellar absorption 
features show the opposite trend. 

Table 2 summarizes various continuum and line 
measurements obtained for the individual regions 
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and for the 3D field of view. We computed the 
H- and if-band flux densities by averaging over 
the spectra of each region in the ranges A = 
1.52 — 1.78 (im and A = 1.96 — 2.42 /zm, respec- 
tively. Throughout this paper, we will adopt the 
photometric system of Wamsteker (1981). 

We measured the emission line fluxes by inte- 
grating over the line profile after subtracting a lin- 
ear continuum fitted to adjacent line- free spectral 
regions. Due to the substantial continuum struc- 
ture, these intervals were carefully selected by in- 
spection of template spectra of K and M stars from 
existing stellar atlases (Klcinmann & Hall 1986; 
Origlia, Moorwood, & Oliva 1993; Dallicr, Bois- 
son, & Joly 1996; Forster Schreiber 2000) so that 
the linear interpolation adequately represented the 
underlying continuum. We tested the validity of 
this procedure by integrating the line fluxes af- 
ter subtracting several template spectra of stars 
within ± 3 spectral classes of the representative 
type for each region (determined in section 5). 
We adopted the average fluxes, with uncertainties 
from the continuum subtraction corresponding to 
one standard deviation of the multiple measure- 
ments. At the spectral resolution of 3D, the Brll 
and Brl2 lines at 1.6807 /xm and 1.6407^m are 
blended with the [Fe II] lines at 1.6769 ^m and 
1.6435 /im, respectively. In these cases, we fitted 
double gaussian profiles to the features. The emis- 
sion lines contribute at most 1.5% to the broad- 
band flux densities for all four regions. 

We measured the equivalent widths (EWs) of 
the Si I absorption feature at 1.59 ^m, and of the 
12 CO(6,3) and 12 CO (2,0) bandheads at 1.62 ^m 
and 2.29 /im (denoted hereafter W1.59, Wi.625 and 
W2.29) according to the definitions and corrections 
for resolution effects given by Forster Schreiber 
(2000). Such corrections were necessary in order 
to compare the EWs for population synthesis pur- 
poses (section 5) with stellar data from existing 
libraries given for R ~ 1600 and R ~ 2500 in the 
H- and K-band, respectively. The Brl4 emission 
line partially fills the Si I feature in all spectra. 
In order to remove this contamination, we scaled 
the profile of the Brl3 line according to the intrin- 
sic ratio Brl4/Brl3 = 0.80 and subtracted it at 
the position of Brl4. This ratio is for case B re- 
combination (Hummer & Storey 1987), assuming 
the electron density and temperature determined 
in section 4 (n c = 300 cm" 3 and T c = 5000 K). 



Differential extinction between these lines can be 
neglected because they are so close in wavelength. 

3.2. Near-infrared Images 

We obtained maps of the broad-band and line 
emission as well as of the EWs from the 3D data 
cubes by applying to each pixel the procedure de- 
scribed above for the spectra. Figure 3 presents 
selected images and ratio maps. 

The spatial distributions of the H- and if-band 
emission obtained with 3D agree well with maps 
in the literature {e.g. Dictz et al. 1986; Tclc- 
sco et al. 1991; Larkin et al. 1994; Satyapal et al. 
1995). The emission is generally centrally concen- 
trated about the nucleus and small-scale structure 
is apparent, in particular the so-called "secondary 
peak" 8" west from the nucleus. We assessed the 
accuracy of the flux calibration and mosaicking by 
comparing the H- and if-band flux densities mea- 
sured at different positions between the 3D data 
and the large-scale images at 1" resolution from 
Satyapal et al. (1997), which were obtained with 
a single detector array. The relative flux densities 
between various regions agree typically to 15% in 
the iJ-band and 10% in the if-band. The differ- 
ences may be due in part to the different spatial 
resolution of the data sets. 

The Br7 and He I 2.058 fim line emission ex- 
hibit clumpy morphologies and follow each other 
very well. The ratio map, however, reveals spatial 
variations in the He I 2.058/Br7 line ratio. The 
contrast between the distribution of the ionized 
gas and broad-band emission is well delineated 
by the logarithmic map of the ratio of if-band 
to Lyman continuum luminosities, log(L k / -^Lyc) 3 • 
The Br7 image obtained with 3D agrees well with 
previously published maps (Waller, Gurwcll, & 
Tamura 1992; Larkin et al. 1994; Satyapal et al. 
1995). Table 3 compares the Br7 line flux inte- 
grated in circular apertures centered on the nu- 
cleus with measurements from the literature. The 
discrepancies between the various results are very 

3 For the if -bandpass A = 1.9 — 2.5 £tm (Wamsteker 1981), 
Lx = 1-87 X 10 19 -D 2 /x where is in Lq (taken as 
3.85 X 10 26 W), D is the distance in Mpc, and f K is the flux 
density in Wm -2 (im -1 . For case B recombination (Hum- 
mer & Storey 1987) with n c = 300 cm" 3 , T c = 5000 K, 
and an average Lyman continuum photon energy of 15 eV, 
Lhyc = 4.68 X 10 22 -D 2 F Br7 where L Lyc is in Lq and the 
Br7 line flux Fbij is in Wm~ 2 . 
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large: up to a factor of 3.5. They are probably 
attributable to uncertainties in absolute calibra- 
tion and continuum subtraction. In addition, posi- 
tioning errors and differences in spatial resolution 
for linemaps could lead to appreciably different 
fluxes given the clumpy morphology of the emis- 
sion. We have however no obvious explanation in 
the case of measurements at higher spectral reso- 
lution from well-registered linemaps, as for Larkin 
et al. (1994) and Satyapal et al. (1995). Both 
the spectral resolution and coverage as well as the 
quality of the 3D spectra allow a better continuum 
subtraction than in previous studies, supporting 
the accuracy of our Br7 fluxes. 

The W2.29 and W1.62 maps reveal very deep ab- 
sorption features around the nucleus and the "sec- 
ondary peak," and progressively shallower features 
at B2 and Bl. The log(Wi.62/W2.29) map is rather 
uniform but indicates an enhancement of W1.62 
relative to W2.29 around Bl. The W2.29 image 
is consistent with the CO index map obtained at 
low spectral resolution by Satyapal et al. (1997) 
within the regions observed with 3D. 

The images demonstrate that the nucleus and 
Bl have the most extreme properties within the 
3D field of view. The 3D maps and spectra in- 
dicate important spatial variations on relatively 
small scales in the composition of the stellar pop- 
ulation, and support that the most recent star 
formation activity — as traced by the Br7 emis- 
sion for example — has taken place outside the 
nucleus, as noted previously (e.g. Telesco et al. 
1991; McLeod et al. 1993; Satyapal et al. 1997). 

3.3. Mid-infrared Spectra 

Figure 4 shows the SWS full scan spectrum 
of M82. Several H recombination lines from the 
Brackett, Pfund, and Humphreys series are de- 
tected as well as pure rotational lines from H 2 
and numerous fine-structure lines from various 
atoms mostly in low ionization stages. Broad 
emission features commonly attributed to poly- 
cyclic aromatic hydrocarbons (PAHs), a broad 
dip centered near 10 /im, and a rising contin- 
uum at A > 10 /im likely due to very small 
dust grains are conspicuous (e.g. Gillett et al. 
1975; Willncr et al. 1977; Lcger, d'Hendecourt, & 
Defourneau 1989; Allamandola, Tielens, & Barker 
1989; Desert, Boulanger, & Puget 1990; Roche et 
al. 1991). According to Sturm et al. (2000), the 



minimum around 10 /im may not be due entirely 
to absorption by interstellar silicate grains as is 
usually assumed; the presence of strong flanking 
PAH emission complexes superposed over a weak 
continuum may contribute to the apparent dip. 
The higher signal-to-noise (S/N) ratio line scans 
are plotted in figure 5. 

We measured the emission line fluxes from both 
the full scan SWS01 spectrum and the individ- 
ual line scan SWS02 spectra. We fitted gaussian 
profiles to the features after subtraction of the 
continuum baseline obtained by fitting a line to 
adjacent line-free portions of the spectrum. Ta- 
ble 4 gives the line fluxes. In order to compare 
the data obtained with the different SWS aper- 
tures, we also scaled the fluxes of the H recombi- 
nation and ionic fine-structure lines to the small- 
est 14" x 20" aperture. We estimated the scaling 
factors from the [Ne II] 12.8 /im map of Achter- 
mann & Lacy (1995); they are 0.8 and 0.7 for the 
14" x 27" and 20" x 33" apertures, respectively, 
with 10% uncertainty. The Bra 4.051 /mi, [Ar III] 
8.99 /mi, and [S IV] 10.51 /mi maps obtained by 
these authors exhibit morphologies similar to that 
of the [Ne II] emission, justifying the use of the 
same scaling factors for all the hydrogen and fine- 
structure lines considered. No beam-size scaling 
was applied to the H2 lines because of the lack of 
information on the spatial distribution of the H 2 
emission. At the resolution of the SWS, the Br/3 
line (2.625 /im) is blended with the H 2 1-0 0(2) 
line at 2.626 /im. However, from the strength of 
the H2 1-0 Q(3) at 2.423 /im, we estimate that the 
H2 1-0 0(2) line contributes at most 30% to the 
measured flux (Black & van Dishoeck 1987; Stern- 
berg & Dalgarno 1989). 

4. NEBULAR ANALYSIS OF M 82 

In this section, we present our nebular analysis 
of M82. We first derive the physical parameters 
critical for the interpretation of the nebular line 
emission: extinction, electron density, gas-phase 
abundances, and ionization parameter appropri- 
ate for the star-forming regions. We then use pho- 
toionization models to constrain the average effec- 
tive temperature of the OB stars. 
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4.1. Interstellar Extinction 

The issue of extinction towards the starburst 
regions of M82 has long been controversial. In 
particular, extinctions in visual magnitudes Ay 
ranging from a few to about 15 mag have been 
obtained under the assumption of a uniform fore- 
ground screen model or from optical and near- 
infrared diagnostics, while Ay w 20 — 60 mag 
have been inferred for a mixed model or from di- 
agnostics at longer wavelengths (see e.g. McLeod 
et al. 1993, Satyapal et al. 1997, and references 
therein). Such differences may substantially af- 
fect the derived intrinsic properties (fluxes, lumi- 
nosities). For instance, the correction factors at 
2.2 /xm range from 1.6 to 4.1 assuming a uniform 
foreground dust screen with Ay = 5 — 15 mag, 
and from 2.2 to 5.7 assuming a mixed model with 
Ay = 20 — 60 mag (see below for the computa- 
tion of the correction factors). The discrepancies 
in extinction estimates can be understood in view 
of the large uncertainties in beam-size corrections 
involved in several studies, of those in the interpre- 
tation of some diagnostics (such as the broad dip 
around 10 ^m discussed above), of the spatially 
non-uniform extinction across the disk of M82, 
and of the large optical depths preventing radi- 
ation at shorter wavelengths to escape from the 
most obscured regions (e.g. Puxley 1991; Telesco 
etal. 1991; McLeod et al. 1993; Larkin et al. 1994; 
Satyapal et al. 1995, 1997; Sturm et al. 2000). 

4-1.1. Global Extinction towards the Ionized Gas 

One major hindrance in previous studies has 
been the lack of data in the mid-infrared regime. 
ISO-SWS observations have now filled this gap, 
with a consistent set including several H recombi- 
nation lines between 3 /xm and 10 /im. H recom- 
bination lines are excellent "standard candles" for 
extinction determinations because their intrinsic 
line emissivities are well determined theoretically. 
Moreover, the derived extinction parameters are 
fairly insensitive to the choice of electron density 
and temperature since the relative line emissivities 
vary slowly with these properties. For instance, 
the results from Hummer & Storey (1987) imply 
that the emissivities of the lines considered below 
change by 5% on average (18% at most) between 
n c = 10 2 cm- 3 and 10 4 cm" 3 , for T c = 5000 K. 
They vary by 20% on average (33% at most) be- 



tween T c = 5000 K and 10 4 K, for n c = 10 2 cnr 3 . 
The intrinsic Lyman continuum photon emission 
rates Q\, yc derived from the dereddened line fluxes 
depend very weakly on n , and vary only slightly 
with T e since the total H recombination coefficient 
«b cx T~ a - 81 . For the examples above, Q^ yc in- 
creases on average by w 5% and 25%, respectively. 

We combined our SWS data with H line mea- 
surements from the radio to the optical regimes, 
obtained mainly in large apertures (diameter 
> 20"), and with appropriate beam-size correc- 
tions. Since the SWS and larger apertures in- 
clude the most prominent emission regions (see the 
linemaps of McCarthy, Heckman, & van Breugel 
1987, Satyapal et al. 1995, Achtermann & Lacy 
1995, and Seaquist et al. 1996), the derived ex- 
tinction should be representative for the bulk of 
ionized gas in M82. A large wavelength cover- 
age is essential for discriminating between various 
dust and sources geometries, since deviations from 
a simple uniform foreground screen model are per- 
ceptible only for diagnostics probing appreciably 
different optical depths. 

In the range A ~ 1 — 10 mm, we used mea- 
surements for H26a, H27a, H41a (Seaquist et 
al. 1996), H30a (Seaquist, Kcrton, & Bell 1994), 
H40a and H53a (Puxley et al. 1989). The in- 
tegrated fluxes for large apertures are consistent 
with predominantly optically thin, spontaneous 
emission in local thermodynamical equilibrium 
(LTE; see references above). We excluded data 
at centimeter wavelengths, which are potentially 
affected by stimulated emission and free-free ab- 
sorption (e.g. Seaquist, Bell, & Bignell 1985; 
Seaquist et al. 1996). At near-infrared wave- 
lengths, we used the Brackctt line fluxes mea- 
sured with 3D. These were complemented with 
the Pa/3 measurements from McLeod et al. (1993) 
and Satyapal et al. (1995), averaging the scaled 
fluxes together because of the large discrepancy 
between the two results (as for Br7 discussed in 
section 3.2). Finally, at optical wavelengths, we 
used the Ha measurement of McCarthy et al. 
(1987). We obtained an additional estimate from 
the Ha + [N II] AA6548, 6584 A flux of Young, 
Kleinmann, & Allen (1988), assuming a uniform 
[N II] /Ha ratio of 0.5 (McCarthy et al. 1987). 

We scaled the line fluxes to match the SWS 
14" x 20" beam in two steps. The fluxes were first 
scaled to a 30"-diameter aperture centered on the 
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nucleus. The beam-size corrections for the mil- 
limeter lines were inferred from those derived by 
Seaquist et al. (1994) between 19", 21", and 41" 
apertures and the entire emission region, with an 
uncertainty of 15%. Those for the Ha measure- 
ments were estimated from the map of McCarthy 
et al. (1987), with 40% uncertainty. For Br7 and 
Pa/3, we applied a scaling oc 9 15 (where 9 is the 
aperture diameter) derived from the Br7 data of 
Satyapal et al. (1995) in 3.8", 8", and 30" aper- 
tures, with 10% uncertainty. Because of the weak- 
ness of the 7/-band Brackett lines within w 5" of 
the nucleus, we used the fluxes integrated over the 
entire 3D map multiplied by the ratio of Br7 fluxes 
scaled to a 30" beam and measured in the 3D 
field of view. In the second step, a scaling factor 
30" -> SWS of 0.5 was applied, with 10% uncer- 
tainty, as derived by comparing various continuum 
and line fluxes from the SWS data with results 
reported in the literature, all of which presum- 
ably trace the same sources since they have similar 
spatial distributions; the references include Klcin- 
mann & Low (1970), Rieke & Low (1972), Gillett 
et al. (1975), Willner et al. (1977), Houck et al. 
(1984), Telesco, Dressel, & Wolstencroft (1993), 
and Achtermann & Lacy (1995). Table 5 lists the 
observed and scaled line fluxes. 

We derived the extinction by least-squares fit- 
ting to 

Fx _( h\( X\ 
F ie{ 



Jref 



»-rcf 



(1) 



where F\/F Te f are the observed line fluxes rela- 
tive to that of a reference line, j\/j Te f are the 
intrinsic ratios of line emissivities, and X\ and 
X rc f are the attenuation factors due to extinction 
at the wavelengths of the lines considered. The 
SWS lines constituting the largest, self-consistent 
data set, we chose Bra as the reference line. We 
took the intrinsic line emissivities from Hummer 
& Storey (1987) for case B recombination with 
n c = 100 cm -3 and T e = 5000 K, appropriate 
for M82 as shown in section 4.2.1 4 . We adopted 
the extinction laws of Draine (1989) and Rieke 
& Lebofsky (1985) at infrared and optical wave- 
lengths, respectively. At millimeter wavelengths, 
the extinction can be neglected. 



4 The electron density is actually in the range ~ 10 — 
500 cm -3 , but we used the tabulated values for 100 cm~ 3 
since the relative line emissivities depend weakly on n c . 



We considered two representative model geome- 
tries. For a uniform foreground screen of dust 
( "UFS" ) , the attenuation is given by 



X UFS 



(2) 



where 7° and 7a are the intensities of the inci- 
dent and emergent radiation, and t\ is the opti- 
cal depth of the obscuring material. For a mixed 
model ( "MIX" ) consisting of a homogeneous mix- 
ture of dust and sources, 



x^MIX _ I* _ 1 e - 

a a - To 



(3) 



Here, 7° is the total intrinsic line intensity pro- 
duced within the mixed medium, and t\ is the 
total optical depth of this medium. The extinc- 
tion in magnitudes is related to the optical depth 
through A\ = 1.086 r A . 

The best fit for each model is achieved with a 
total visual extinction of A V FS = A+\' 5 mag and 
Ay lx — 43 ± 23 mag. We estimated the uncer- 
tainties on Ay from the reduced chi-squared dia- 
grams, with a la error corresponding to a factor 
of e 1 / 2 from the minimum x\- Table 6 compares 
the best-fit extinction-corrected line fluxes rela- 
tive to Bra for each geometry (columns labeled 
"Draine") with the intrinsic line ratios. Figure 6 
illustrates the results (plots labeled "Draine"), 
as the Ql vc derived from each of the extinction- 
corrected line fluxes assuming an average Lyman 
continuum photon energy of 15 eV. The average 
values of Ql yc are given in table 6. 

Good extinction models are those for which 
equal values of Q^yc are inferred from each of the 
extinction-corrected line fluxes. Figure 6 demon- 
strates that purely foreground extinction provides 
a much less satisfactory fit to the data over the 
entire wavelength range considered, as also shown 
by e.g. Puxlcy (1991) and McLeod et al. (1993). 
The best fit for this geometry results in a minimum 
Xn ~ 8 compared to x n ~ 2 for the mixed model. 
For the latter geometry, the extinction corrections 
are of 4.8 near 2 /zm and 1.6 near 5 fim, with un- 
certainties of w 50% and « 20%, respectively. 

Increasing n c to 10 4 cm~ 3 or T e to 10 4 K re- 
sults in larger Ay's, but within the la uncertain- 
ties for the nominal case. The implied values of 
Ql yc increase by about 15% and 50%, respectively. 
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We also performed fits excluding the Ha measure- 
ments, which may include a component from light 
escaping along the minor axis scattered by dust 
grains (e.g. O'Connell & Mangano 1978; Notni 
1985). The Ha flux could thus be overestimated 
relative to the data at longer wavelengths, since 
the scattering efficiency of interstellar dust grains 
generally decreases rapidly with increasing A (e.g. 
Emerson 1988). The results for the mixed model 
are little affected, with Ay IX = 44 mag (a sub- 
stantially larger A V FS = 9 mag is obtained, as 
expected since the dust and sources are actually 
mixed). Fits excluding the millimeter lines yield 
A™ IX = 18 mag, illustrating the importance of 
including unobscured lines. 

4-1-2. The A = 3 — 10 /im Extinction Law 

Despite the good overall fit for the mixed 
model, figure 6 indicates some deviations in the 
A = 3 — 10 /im region. The constraints imposed by 
the millimeter lines prevent reduction of these de- 
viations by a change in absolute level of extinction. 
The extinction law assumed for these wavelengths 
may not be appropriate however. Indeed, until 
recently, the 3 — 10 /im extinction law was poorly 
determined because of the difficulties inherent to 
ground-based observations and because the prop- 
erties of the template sources accessible so far in 
this range were not well-known. The SWS has 
now provided observations of numerous nebular H 
recombination lines between 3 [im. and 10 ^m in a 
variety of objects. These lines have been used to 
investigate the extinction law in the direction of 
the Galactic Center (Lutz 1999). The "Galactic 
Center extinction law" (hereafter simply GC law) 
lacks the pronounced minimum in the 4 — 8 ^tm 
region expected for standard graphite-silicate dust 
mixtures which are usually assumed (e.g. Draine 
1989 and references therein), suggesting additional 
contributors to the extinction. 

The SWS line fluxes in M82 arc much bet- 
ter reproduced if the GC law is adopted. The 
best fits are obtained with ^4y FS = 4+3 mag and 
Ay lx — 52 ± 17 mag. The results are given in 
table 6 and figure 6 (columns and plots labeled 
"GC"). The fit for the uniform foreground screen 
model is still much poorer than for the mixed ge- 
ometry (xi ~ 8 and « 1, respectively). The ex- 
tinction corrections for the mixed model are 5.8 
near 2xtm and 2.6 near 5xtm, with uncertainties 



of « 30% and 25%, respectively. Although the 
data do not allow the accurate determination of 
the extinction law in M82, they provide evidence 
for deviations from the commonly used Draine law 
similar to those found towards the Galactic Center 
given by Lutz (1999). 

4-1.3. Adopted Parameters 

In view of the above analysis, we will adopt 
the mixed model with Ay — 52 mag through- 
out this paper as representative of the global ex- 
tinction towards the bulk of the ionized gas in 
M82. The Lyman continuum photon emission 
rate for a 30"-diameter aperture is twice that for 
the SWS field of view (section 4.1.1), Q° yc = 
1.23 x 10 54 s -1 . Our results are consistent with 
the Q^yc inferred from the millimeter thermal free- 
free emission (~ 10 54 s _1 ; e.g. Carlstrom & Kro- 
nberg 1991; Seaquist et al. 1996), and confirm 
the results of McLeod et al. (1993) who found 
A™ IX = 55 mag and Q£ yc = 1.05 x 10 54 s" 1 . 
We will also adopt the GC extinction law of 
Lutz (1999) between 3 /zm and 10 /xm, keeping the 
Draine (1989) law for the other relevant infrared 
ranges and the Ricke & Lebofsky (1985) law at op- 
tical wavelengths. Satyapal et al. (1995) demon- 
strated the validity for M82 of the Draine law 
at near-infrared wavelengths (A\ oc A -175 ), or of 
similar extinction laws (e.g. Landini et al. 1984, 
with A x oc A" 1 - 85 ). 

4-1-4- Local Extinction towards the Ionized Gas 

From the spatially resolved Brackett line emis- 
sion obtained with 3D, we derived the extinction 
towards individual regions as described in sec- 
tion 4.1.1 using Br7 as the reference line. The 
results for selected regions are reported in table 7, 
together with the <3L yc computed from the Br7 
fluxes (the most accurately measured line), cor- 
rected for Ay lx . Due to the weakness of the H- 
band lines over significant areas, we could only 
generate a partial extinction map, from the 3D 
lincmaps rebinncd to 1" x 1" pixels to increase 
the S/N ratio. 

The variations in Ay between individual re- 
gions and in our partial extinction map agree well 
with those seen in the extinction map of Satyapal 
et al. (1995) obtained from Br7 and Pa/3 measure- 
ments. The lower Ay IX for the 3D field of view 
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compared to the global extinction derived above 
probably reflects the optical depth limitations of 
the near-infrared diagnostics. For the same rea- 
sons, the derived Ay IIX and Ay FS imply similar 
extinction corrections; for example, the differences 
are 35% or less near 2 /mi for the selected regions, 
and < 20% on the scale of the rebinncd pixels over 
the valid regions. The uncertainties on the extinc- 
tion corrections are typically 20% — 40% for both 
models. As will be discussed in section 4.3, the 
ionized nebulae arc likely mixed with the molecu- 
lar gas and dust clouds even on scales of a few tens 
of parsecs. We will therefore adopt the results for 
the mixed model. 

4.2. Physical Conditions of the ISM 

4-2.1. Electron Temperature and Density 

From the ratios of H recombination lines with 
the (thermal) continuum flux density in the mil- 
limeter regime, a T c near 5000 K has been derived 
for the starburst core of M82 (Puxlcy et al. 1989; 
Carlstrom & Kronberg 1991; Seaquist et al. 1994, 
1996). From the results of Seaquist et al. (1996), 
based on maps of the H41a line emission and of the 
underlying continuum, we inferred that this tem- 
perature is appropriate on smaller spatial scales 
for most of the regions observed with 3D. 

The SWS spectrum provides three density- 
sensitive ratios: [S III] 18.7/mi/33.5 /mi, [Ne III] 
15.6 ^m/36.0 /mi, and [Ar III] 8.99/mi/21.8/mi. 
All three are fairly insensitive to T e in the range 
5000 - 20000 K. The [S III] ratio is the most 
reliable one because both lines are amongst the 
strongest in the SWS spectrum. It is also the 
most sensitive at low n c because the upper levels 
of the transitions have the lowest critical den- 
sities. For the 14" x 20" SWS aperture, the 
dereddened suplhur, neon, and argon ratios are 
0.71 ± 0.23, 9.0 ± 4.2, and 15.5 ± 9.3, respectively. 
The uncertainties include those of the line flux 
measurements (from the relative flux calibration 
and continuum subtraction), and of the beam- 
size and extinction corrections. Comparison with 
results of computations of collisional excitation 
(e.g. Alexander et al. 1999) shows that the mea- 
sured ratios lie in the low-density limit, indicating 
n e ~ 10 - 600 cm- 3 . 

Our result is similar to the average n c in large 
apertures obtained previously by various authors 



using different infrared and radio diagnostics, and 
assuming a single-density model (Houck et al. 
1984; Duffy et al. 1987; Seaquist et al. 1985, 
1996; Colbert et al. 1999). We also considered 
the [O III] 52/mi/88/mi ratio, which is more sen- 
sitive to n c at low densities. The ratio measured 
by Duffy et al. (1987) in a 48"-diameter aperture, 
similar to the value for an 80"- aperture of Col- 
bert et al. (1999), indicates n c <~ 50 — 500 cm~ 3 . 
We will adopt a value of 300 cm~ 3 . 

4-2.2. Gas-phase Abundances 

We determined the gas-phase abundances of 
Ne, Ar, and S, three of the most abundant heavy 
elements in H II regions, using the lines detected 
with the SWS. These have critical densities for col- 
lisions with electrons at T = 5000 K in the range 
~ 10 — 10 6 cm" 3 so that collisional de-excitation 
can be neglected. Assuming a "one-layer" model 
with uniform density and temperature, the ionic 
abundances can be computed from 

F A(X+») = "(X+-)"cjA(X+-) ^ 
F\(H+) ™(H+)«c/A(H+) 

where i*A(x+*) and -F\(h+) are the fluxes of the 
ionic line of interest and of a reference H recom- 
bination line, n(x+i) and n(H+) are the densities 
of ions X +1 and H + , and jA(x+ i ) and Ja(h+) are 
the line emissivities. In H II regions, H is nearly 
completely ionized so that n H + s=a njj. We com- 
puted the fine-structure line emissivities in the 
low-density limit using the effective collisional 
strengths from Johnson, Kingston, & Dufton 
(1986), Saraph & Tully (1994), Butler & Zeippen 
(1994), Pelan & Berrington (1995), and Galavi's, 
Mendoza, & Zeippen (1995). We took Bra as ref- 
erence H line, with its emissivity from Hummer & 
Storey (1987). 

Table 8 gives the data and results. For ioniz- 
ing stars with effective temperatures in the range 
35000 K - 40000 K, as found in section 4.4 below, 
the Ne, Ar, and S are expected to be mostly in 
the ionization stages observed with SWS. The ele- 
mental abundances are thus well approximated by 
the sum of the ionic abundances determined here. 
The abundances are nearly solar or slightly above 
for Ne and Ar, and about one-fourth solar for S. 
Similar underabundances for S have been found in 
Galactic H II regions (e.g. Simpson et al. 1995) 
and in some extragalactic starburst systems (e.g. 
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Genzel et al. 1998), and are attributed to deple- 
tion of S onto interstellar dust grains. 

Most abundance determinations for M82 in 
the literature indicate no large depletions or en- 
hancements for most elements compared to the 
solar neighbourhood composition (e.g. Gaffney & 
Lester 1992; McLcod et al. 1993 and references 
therein; Achtermann & Lacy 1995; Lord et al. 
1996). The exception is Si, for which Lord et al. 
found a gas-phase abundance three times larger 
than in Galactic nebulae, and which they interpret 
as probably due to partial destruction of silicate 
grains by fast supernova-driven shocks. 

4.3. Ionization Parameter 

For the purpose of photoionization modeling, 
we represented H II regions as thin gas shells sur- 
rounding central, point-like stellar clusters (see 
top of figure 7). In such "central cluster" mod- 
els, the nebular conditions are specified by the 
distance R between the ionizing cluster and the 
illuminated surface of the gas shell, the H num- 
ber density nn, and the ionization parameter U 
defined as 

jj _ Qhyc _ 0Lyc 

where QLyc is the production rate of Lyman con- 
tinuum photons from the stars and c is the speed 
of light. U thus gives the number of Lyman con- 
tinuum photons impinging at the surface of the 
nebula per H atom. Since H is nearly completely 
ionized in H II regions, and since He is not fully 
ionized in M 82 (section 4.4), we assumed tir ~ n e . 

In such complex and distant systems as star- 
burst galaxies, a large number of H II regions may 
coexist in a relatively small volume and may not be 
individually resolved by the observations. Thus, 
in constraining U from the observed properties, 
the shell geometry may not be directly applica- 
ble. The spatial distribution of the gas relative to 
the sources is crucial in determining the Lyman 
continuum photon flux </>L yc impinging on the gas, 
and in deriving the effective U (hereafter U e s) to 
model appropriately the nebular emission in the 
framework of the idealized central cluster geom- 
etry. Due to the lack of information on small 
enough spatial scales, the ionization parameter is 
generally poorly determined in starburst galaxies. 
M 82 is one exception owing to its proximity. We 



have used our 3D and SWS data together with 
data from the literature to constrain the degree of 
ionization of the nebulae within the starburst re- 
gions. Appendix A gives the details of our deriva- 
tion; in the following, we restrict ourselves to out- 
lining the main results. 

The ISM properties within the starburst regions 
of M82 (e.g. Lord et al. 1996) suggest it can be 
represented by a collection of clouds with, on aver- 
age, a molecular core with radius ~ 0.5 pc shielded 
by a thin neutral atomic layer, and an ionized layer 
extending out to « 1 pc. These clouds have a 
mean separation of « 2 — 7 pc. High resolution op- 
tical and near-infrared imaging reveals large num- 
bers of young, luminous clusters throughout the 
starburst core (e.g. O'Connell et al. 1995; Satya- 
pal et al. 1997). Assuming the OB stars reside in 
such clusters and adopting a plausible cluster lu- 
minosity function consistent with the properties of 
the optically-selected population, the intrinsic Ly- 
man continuum photon emission rates imply an 
average cluster separation of 2 — 9 pc. 

The similar cloud-cloud and cluster-cluster 
mean separations suggests that the clouds and 
clusters are well-mixed and uniformly distributed 
throughout the regions of interest, as illustrated 
at the bottom of figure 7. For such a random 
distribution model, the flux of Lyman continuum 
photons impinging on the nebulae is reduced com- 
pared to the central cluster geometry due to the in- 
crease in surface area of the gas exposed to the ra- 
diation field. From the analysis of selected regions 
in appendix A, we determine log [/eg- w —2.3 dex 
and conclude that it is representative of the lo- 
cal nebular conditions throughout all of the star- 
forming regions in the central 500 pc of M 82. We 
also show that U e s can be expressed in terms of 
the ratio of intrinsic Lyman continuum photon 
emission rate and of the molecular gas mass, with 
U c ff oc Qhyc/M}i 2 and the proportionality factor 
depending only on the gas cloud properties (see 
equation [A4]). Since the cloud properties are 
similar throughout the starburst core of M82, the 
comparable values of log U c s are interpreted as re- 
sulting from comparable star formation efficiencies 
as measured by Q-Lyc/Mn 2 . 

4.4. Young Stellar Populations 

The SWS and 3D data provide several diagnos- 
tics sensitive to the shape of the ionizing radia- 
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tion spectrum, dominated by OB stars in M82. 
These include ratios of mid-infrared atomic fine- 
structure lines and of near-infrared He to H recom- 
bination lines, which probe the 13 eV — 41 eV en- 
ergy range. The contribution from shock-ionized 
material to the line emission considered below is 
not likely to be important in M 82 (McLeod et al. 
1993; Lutz et al. 1998). In the following, we will 
assume that the lines originate entirely in gas pho- 
toionized by the OB stars. Table 9 summarizes the 
data and results. 

4-4-1- Mid-infrared Fine-structure Line Ratios 

Amongst the diagnostic line ratios available 
from the SWS data set, we considered [Ne III] 
15.6^m/[Ne II] 12.8 /m, [Ar III] 8.99^m/[Ar II] 
6.99 /xm, and [S IV] 10.5/zm/[S III] 18.7,uni to 
avoid complications due to the uncertainties of the 
elemental abundances. 

We modeled the variations of the line ratios 
with effective temperature of the stars (T^ B ) using 
the photoionization code CLOUDY version C90.05 
(Fcrland 1996) and the stellar atmosphere models 
for solar-metallicity main-sequence stars of Paul- 
drach et al. (1998). Wc adopted the nebular pa- 
rameters derived in the previous subsections: solar 
gas-phase abundances, n e = 300 cm~ 3 , and an ef- 
fective log U — — 2.3 dex. With this value for the 
ionization parameter, equation (5) implies a corre- 
sponding R of several tens to several hundreds of 
parsecs, depending on the region. Since photoion- 
ization models are little sensitive to variations in R 
above ~ 10 pc, we adopted a fixed R = 25 pc. We 
neglected the effects of dust grains possibly present 
within the nebulae. Dust grains mixed with the 
ionized gas are not expected to affect significantly 
the ionization equilibrium of H II regions because 
the variations of the dust absorption cross-section 
with wavelength resembles closely that of H, peak- 
ing near 17 eV (e.g. Mathis 1986). Mathis further 
argues that this conclusion is not altered for dust 
properties and gas-to-dust ratios characteristic of 
the Galactic and Magellanic Clouds diffuse ISM or 
of denser H II regions. 

Figure 8 shows the theoretical predictions and 
the measured ratios, corrected for extinction and 
beam-size differences when appropriate. The ef- 
fects of variations of n e and log U, or of adopting a 
gas and dust composition typical of the Orion neb- 
ula are also indicated. The most sensitive parame- 



ter affecting the line ratios is the ionization param- 
eter. However, varying log U in the plausible range 
from —2 to —2.5 dex (see appendix A) implies 
relatively small differences in Tjg 8 : < ±1000 K 
for the neon and argon ratios, and —2000 K or 
+4000 K for the sulphur ratio. Given the uncer- 
tainties of the data (line measurements, extinc- 
tion and aperture corrections) as well as of the 
models (e.g. nebular parameters, stellar atmo- 
spheres, atomic data), the agreement between the 
results from the three diagnostic ratios is satis- 
factory. Our high quality and consistent data set 
confirms results obtained in the past from various 
infrared and millimeter lines, which ranged from 
30000 K to 37000 K (Gillett et al. 1975; Willner et 
al. 1977: Puxlcy et al. 1989; McLeod et al. 1993; 
Achtermann & Lacy 1995; Colbert et al. 1999). 

4-4-%- Near-infrared He to H Recombination Line 
Ratios 

The strongest recombination lines from singly- 
ionized He detected in the 3D spectra correspond 
to the 2 1 S" - 2 1 P triplet transition at 2.058 /jm 
and the 3 3 P — 4 3 D singlet transition at 1.701 /im 
5 . Together with the nearby Br7 and BrlO lines, 
respectively, they provide T^ B diagnostics sensi- 
tive in the range < 55000 K and very little affected 
by extinction. The low-excitation SWS spectrum 
and the non-detection of He II lines in the 3D data 
(for example at 2.189 ^m) rule out the presence 
of important populations of Wolf-Rayet stars or 
other sources that are hot enough to doubly ion- 
ize He, which would complicate the interpretation 
of the ratios. 

The He I 2.058/Br7 line ratio has been com- 
monly used to estimate T e g B in near-infrared stud- 
ies of starburst galaxies and Galactic H II regions 
(e.g. Doyon, Puxley, & Joseph 1992; Doherty 
et al. 1994, 1995). However, due to sensitivity 
to local physical conditions and to degeneracy in 
TtfP through resonance and collisional effects (e.g. 
Robbins 1968; Clegg 1987; Shields 1993), this ra- 
tio alone is not sufficient for constraining T C ^ B . On 
the other hand, the 1.701 fjxa line originates from a 
higher quantum state in the n 3 P—n' 3 D series and 
is essentially unaffected by self-absorption and col- 
lisional effects but the He I 1.701/BrlO ratio satu- 

5 Wc thank M. and G. Rickc for drawing our attention to the 
He I 1.701 fim line. 



12 



rates above T^P w 40000 K (see references above). 
Therefore, the combination of the He I 2.058/Br7 
and He I 1.701/BrlO ratios allows one to discrimi- 
nate between the high- and low-T^ 6 regimes (see 
e.g. Vanzi et al. 1996 and Doherty et al. 1995 for 
earlier applications). 

We modeled the He to H line ratios using 
CLOUDY and the same model atmospheres and 
nebular parameters as for the mid-infrared line 
ratios. CLOUDY computes the 2.058 /jm line 
but not the 1.701 yum line which we derived in- 
directly as explained below. While models of 
He I 2.058/Br7 presented by various authors in the 
past generally agree very well in the range T e g < 
40000 K, the results at higher T C ° B vary impor- 
tantly (Doyon et al. 1992; Shields 1993; Doherty 
et al. 1994, 1995; Langon & Rocca-Volmerange 
1994). In this regime, the ratio is particularly sen- 
sitive to the treatment of the He I Lya opacity 
which is one of the major sources of uncertain- 
ties in theoretical predictions. Ferland (1999) dis- 
cusses the revised treatment implemented in the 
CLOUDY version we have used. 

For He I 1.701 ^m, we used the He I 4471 A flux 
predicted by CLOUDY. Both lines originate from 
the 4 3 D level and their fluxes are proportional 
to each other since triplet transitions from high 
n 3 D levels are essentially unaffected by collisional 
effects or scattering from the metastable 2 3 S level. 
We scaled the resulting He I 447lA/BrlO curve 
so that its saturation value for full He ionization 
equals that for He I 1.701/BrlO given by Vanzi et 
al. (1996). While n c is fixed in running CLOUDY, 
T c varies for each T e g- B ; we thus interpolated the 
relationship of Vanzi et al. in T c as appropriate. 
Since collisions are unimportant for both He lines 
involved here, we used this relationship up to n c = 
10 3 cm" 3 . 

The ratios measured for selected regions are 
plotted against the model predictions in figure 8. 
The He I 1.701/BrlO ratio clearly rules out the 
high-T c ^ B solutions from the He I 2.058/Br7 ra- 
tio. For the temperatures near 36000 K inferred, 
He I 2.058/Br 7 and He I 1.701/BrlO arc little af- 
fected by variations in n e or log U within plausible 
ranges for M82, or by modest changes in the gas 
and dust composition. Again, in view of the un- 
certainties of the data and models (in particular 
the continuum subtraction for the -ff-band lines), 
the T^'s inferred for each region from both ratios 



agree satisfactorily. 

The variations in our He I 2.058/Br7 map (fig- 
ure 3) support a general though small increase 
in T^P from the nucleus to larger projected radii 
along the galactic plane of M 82 to the west. Such 
radial variations have been suggested by Satyapal 
et al. (1995) from the ratio of Br 7 to 3.29 (im PAH 
feature emission. McLeod et al. (1993) also sug- 
gested such a trend from the lower temperatures 
derived from mid-infrared diagnostics compared to 
those inferred from optical diagnostics assuming 
the former trace the innermost stellar population 
while the latter, foreground and thus outermost 
clusters. On the other hand, Achtermann & Lacy 
(1995) concluded that T^ B is roughly constant 
across the starburst core from the lack of clear 
spatial variations in the excitation state from their 
Bra, [Ne II] 12.8 ^m, [Ar III] 8.99 jum, and [S IV] 
10.5 /zm maps. 

We constrained quantitatively the spatial vari- 
ations in T^p across the 3D field of view from 
our He I 2.058 /xm and Br7 linemaps, rebinned to 
1" x 1" pixels. We neglected extinction effects, 
but this introduces errors < 5% in the ratios. The 
resulting T c % B 's vary from 33200 K to 37500 K, 
with an average of 35700 K and a relatively small 
dispersion of la = 650 K. The 3D data indicate 
therefore a roughly constant T g B for the hot mas- 
sive stars across the regions observed, with only a 
marginal gradient with projected radius. 

Jf..Jf..3. Additional Remarks 

From our data sets, the neon ratio is probably 
the most reliable indicator for the absolute T q B . 
Both lines involved are strong in M 82 and, among 
the three mid-infrared line pairs, their proximity 
in wavelength minimizes most extinction effects 
and they were observed through the same aper- 
ture. It is also the most sensitive of all our diag- 
nostics, probing the largest range in ionizing en- 
ergy (between 22 eV and 41 eV). In addition, the 
inferred T e g B is little affected by the uncertain- 
ties on the physical conditions within the nebu- 
lae. The He I 2.058/Br 7 and He I 1.701 /Br 10 
are potentially less reliable indicators for the ab- 
solute T^P because directly sensitive to the He 
abundance. On the other hand, for the ranges ob- 
served in M82, these ratios are little affected by 
modeling uncertainties and are thus robust diag- 
nostics for the relative variations in T e g B (assum- 
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ing negligible gradients in He abundance). Since 
the smaller SWS aperture and the 3D field of view 
cover the same most prominent sources, the neb- 
ular line emission from both data sets traces es- 
sentially the same stellar populations. In the rest 
of this work, we will adopt the result from the 
neon ratio (37400 K) for the 3D and SWS fields of 
view and will apply a correction of +1400 K to the 
temperatures inferred from He I 2.058/Br7, corre- 
sponding to the difference in obtained from 
the neon ratio and that for the 3D field of view. 

From the temperature scale of Vacca, Garmany, 
& Shull (1996), the inferred T c % B 's within the 3D 
and SWS fields of view correspond to 08.5 V stars, 
with Q£ y 8 c 5V = 10 48 - 72 s" 1 . Given the roughly 
constant nebular excitation and the SWS aper- 
ture including half of the integrated emission in 
a 30"-diameter region centered on the nucleus of 
M82, this spectral type is probably representative 
of the dominant OB stars for the starburst core as 
well. The number of equivalent 08.5 V stars re- 
quired to produce the intrinsic Lyman continuum 
luminosity in various regions is given in table 9. 

5. STELLAR POPULATION SYNTHE- 
SIS OF M 82 

In this section, we apply population synthesis to 
our 3D data. We constrain the spectral type and 
luminosity class of the evolved stars, and investi- 
gate their metallicity. We also constrain the con- 
tribution from additional continuum sources (hot 
dust, OB stars, nebular free- free and free-bound 
processes) as well as the extinction towards the 
evolved stars. All results are reported in table 10. 

5.1. Analysis of Selected Absorption Fea- 
tures 

5.1.1. The Giants- supergiants Controversy 

The nature of the evolved stellar population in 
M82 has long been debated. Various diagnostics 
have been used in the past, including the near- 
infrared broad-band colours, the CO p h and -f^Oph 
photometric indices measuring the depth of the 
CO bandheads longwards of 2.3 /<m and of the 
H 2 absorption feature at 1.9 /zm, spectral syn- 
thesis in the range 2.18 /im — 2.28 /im, and mea- 
surements of the ratio of stellar mass to intrinsic 
if-band luminosity M*/Lk (Walker, Lebofsky, & 
Rieke 1988; Lester et al. 1990; Gaffney & Lester 



1992; Gaffney, Lester, & Telesco 1993; McLeod et 
al. 1993; Lancon, Rocca-Volmerange, & Thuan 
1996). Evolutionary synthesis models have also 
been applied to obtain indirect constraints (Rieke 
et al. 1980, 1993; Satyapal et al. 1997). 

However, no consensus has been reached yet, 
in particular concerning the nucleus: some of the 
above studies indicate the presence of young su- 
pergiants while others provide evidence for old, 
metal-rich giants as dominant sources of the near- 
infrared continuum emission. Possible causes for 
these discrepancies include positioning uncertain- 
ties and dependence on aperture size due to impor- 
tant spatial variations in the indicators, difficulties 
inherent to measurements of H 2 O p h (in a range of 
poor atmospheric transmission), and the weakness 
of several absorption features. More importantly, 
the diagnostics used so far exhibit degeneracy in 
temperature and luminosity, and are affected by 
extinction and featureless continuum emission. 

The 3D data allow us to apply alternative diag- 
nostic tools. The CO bandheads at 2.29 /im and 
1.62 /im together with the Si I feature at 1.59 /im 
are particularly useful in stellar population stud- 
ies from moderate-resolution near-infrared spec- 
tra (Origlia et al. 1993; Oliva et al. 1995; Forster 
Schrcibcr 2000). They suffer much less from mea- 
surements uncertainties and provide sensitive indi- 
cators for the effective temperature and luminosity 
class of cool stars. Moreover, their EWs (W2.29, 
W1.62, and VF1.59) are independent of extinction 
and provide a means of constraining the contri- 
bution, or "dilution," from featureless continuum 
emission sources without requiring any assump- 
tions on their nature and physical properties. 

In the following analysis, we neglect possible 
contributions from thermally-pulsing AGB stars 
(such as Mira variables and N-type carbon stars) 
since none of their extreme, characteristic features 
(e.g. Johnson & Mendez 1970; Langon et al. 
1999) are seen in the 3D spectra. Moreover, evo- 
lutionary synthesis models show that these stars 
never produce more than 10% — 40% of the inte- 
grated near-infrared light of stellar clusters of any 
age (e.g. Bruzual & Chariot 1993; paper 2). 

5.1.2. Selected Regions 

We analyzed the EWs using the diagnostic dia- 
grams proposed by Origlia et al. (1993) and Oliva 
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et al. (1995). These are shown in figures 9 and 
10, where the stellar data has been obtained from 
existing relevant libraries (as compiled by Forster 
Schreiber 2000). The horizontal bars indicate the 
measurements for selected regions in M82; those 
for the 3D field of view are omitted for clarity, but 
are very close to those for B2. Figure 9 gives the ef- 
fective temperature (T e g) and luminosity class im- 
plied by the EWs. Figure 10 allows the determina- 
tion of the amount of dilution near 1.6 /an (D\,q) 
from the vertical displacement relative to the locus 
of stars in the W\.Q2 versus log(Wi. 62/W1. 59) dia- 
gram, and near 2.3 /im (-D2.3) from the horizontal 
displacement in the \V1.e2 versus log(M / i.62/W / 2.29) 
diagram once W1.Q2 is corrected for dilution. As 
shown by Oliva et al. (1995), undiluted composite 
stellar populations fall on the distributions defined 
by the stars in these plots. 

For the central 35 pc at the nucleus of M82, 
both W1.62 and log(Wi.62/Wi.59) indicate an av- 
erage effective temperature in the range 3600 — 
4100 K, implying negligible dilution around 
1.6 /im. The VF2.29 is characteristic of either gi- 
ants with T c ff = 3200 — 3600 K or supergiants 
with T cff = 3600 - 4200 K. Therefore, the EWs 
can only be reconciled for a population of su- 
pergiants and negligible dilution in both H- and 
if -band. Using the temperature calibration of 
Schmidt-Kalcr (1982), the corresponding average 
spectral type is K5 I. For the other selected re- 
gions, degeneracy in luminosity class and dilution 
complicates the interpretation of the EWs. In all 
cases, Wi.62 and log(Wi.62/Wi.59) imply negligible 
dilution near 1.6 /im. For Bl, the data indicate 
dominant K3 III stars and D 2 .z ~ 10% or K2 I 
stars and D 2 . 3 ~ 10% - 40%. For B2 and the 3D 
field of view, the EWs are more consistent with 
K3-K4 supergiants and small amounts of dilution 
near 2.3 /im, but undiluted emission from K4 gi- 
ants cannot be completely ruled out. However, 
we favoured the supergiants solutions on the basis 
of the ratio of stellar mass to intrinsic if -band 
luminosity M* /Lk, as explained below. 

5.1.3. Spatially Detailed Analysis 

We applied a similar analysis to all the regions 
observed with 3D, using maps of the VF1.59, W\.q 2i 
and W2.29 generated from the data cubes rebinned 
to 1" x 1" pixels. We corrected the W1.59 map for 
contamination by the Brl4 emission line as de- 



scribed in section 3.1. The resulting spectroscopic 
indices for all pixels are plotted in the diagnostic 
diagrams of figure 10. 

Together with the EW maps from figure 3, fig- 
ure 10 reveals variations in the intrinsic composi- 
tion of the evolved stellar population on small spa- 
tial scales, little dilution around 1.6 /zm, and vari- 
able dilution around 2.3 /im. The inferred T e fj's 
range from 4500 K down to 3600 K, corresponding 
to spectral types G9 to M0 for supergiants. The 
average is 4000 K with dispersion of la = 200 K, 
equivalent to K4 ± two spectral sub-classes assum- 
ing supergiants. Several regions lie on the locus of 
supergiants in the W1.62 versus log(Wi. 62/1^2. 29) 
diagram; they also coincide with the brightest K- 
band sources. Others are characterized by too 
small 1^2.29 relative to W\.q 2 compared to normal 
evolved stars, implying significant dilution near 
2.3 /im up to w 50% assuming an intrinsic pop- 
ulation of supergiants (« 20% for giants). Be- 
cause dilution in the ii-band is negligible, the 
W1.62 ma P constitutes essentially a T e s map for 
the evolved stars, showing spatial variations that 
are more complex than simple radial gradients. 
The coolest populations are found around the nu- 
cleus and along a ridge extending up to the sec- 
ondary if -band peak (w 8" to the west). Just 
south from this ridge, the T e g increases progres- 
sively along the Nucleus — > B2 — > Bl sequence. 

For most individual regions, the log {W\. q 2 / W 2 . 29) 
does not allow the discrimination between giants 
and supergiants due to degeneracy in luminos- 
ity class and dilution. They correspond predom- 
inantly to the smoother low-surface brightness 
regions in the 3D broad-band maps. The anal- 
ysis of the M* /Lk ratio provides an additional 
constraint. From the stellar mass derived in 
appendix B, the M*/Lk ratio within the star- 
burst core of M82 is very low (« 1.4 M Q /L Q 
and w 0.5 M©/L Q for the central 35 pc and 
500 pc, respectively) 6 . From the 3D maps and 
the data of Satyapal et al. (1997), we estimate 
that the faint smooth emission component repre- 
sents about 75% of the total intrinsic Lk within 
the central 30". This implies an upper limit of 
M* /Lk ~ 0.6 M Q /L Q for the corresponding pop- 

6 We remind the reader that we are using the definition 
L K [L Q ] = 1.87 x 10 19 (D [Mpc]) 2 (f K [Wm- 2 H) for 
the iiT-band covering A = 1.9 — 2.5 /j,m in the photometric 
system of Wamsteker (1981), and with L = 3.85 X 10 26 W. 
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ulation assuming it contains all the mass. 

The above ratios are substantially lower than 
found for old populations in elliptical galaxies and 
bulges of spiral galaxies ("normal populations"), 
which lie typically in the range 10 — 30 M Q /L 
(e.g. Devereux, Becklin, & Scoville 1987; Oliva 
et al. 1995; Hunt et al. 1999). Alternatively, 
red giants, which have a characteristic M/Lk ~ 
0.1 M Q /L , would contribute <~ 15% of the to- 
tal mass if they dominate the low-surface bright- 
ness emission. This is inconsistent with the typical 
fraction of 0.2% — 1% determined empirically for 
normal populations (e.g. Pickles 1985). These 
arguments strongly suggest that young red su- 
pergiants dominate the near-infrared continuum 
throughout the entire starburst core of M82. 

The important and complex spatial variations 
in the composition of the evolved stellar popula- 
tion revealed by 3D indicate that studies based on 
data obtained through different apertures or at a 
few positions only may be misleading. In particu- 
lar, from a comparison of the CO p h measured at 
the nucleus and at the secondary if-band peak, 
Lester et al. (1990) and McLeod et al. (1993) 
concluded that there are no gradients in the com- 
position of the stellar population across the star- 
burst core of M82. However, the 3D EW maps 
clearly show that these regions are "privileged" in 
the sense that they sample populations with very 
similar properties. 

From the EWs at the nucleus, a dominant pop- 
ulation of red supcrgiants is inferred down to the 
central « 20 pc. The spatial resolution of the 3D 
images prevents a reliable investigation of the nu- 
clear population on smaller scales. Gaffney et al. 
(1993) compared their M*/L K ratio in M82 to 
that in the Galactic Center, both within the same 
radius of 7.5 pc, and argued that their similar- 
ity supports an old bulge population as dominant 
nuclear if -band source in M82. We simply note 
here that the near-infrared light within 7.5 pc of 
the Galactic Center contains a significant contri- 
bution from red supergiants (e.g. Hallcr & Rieke 
1989; Blum, Scllgren, & DePoy 1996). 

5.1.4- The Metallicity of the Evolved Stars 

The analysis of the EWs presented above 
is based on empirical indicators valid for stars 
with near-solar metallicities. This seems justified 



given the gas-phase abundances derived in sec- 
tion 4.2.2. The metallicity of the evolved stars 
can in fact be directly constrained using the di- 
agnostics proposed by Origlia et al. (1997) and 
Oliva & Origlia (1998). These are based on Wi.62 
and IL2.29, and derived from theoretical model- 
ing of their behaviours with T e fj, surface gravity, 
micro-turbulent velocity, and metallicity. 

We first examine the central 35 pc of M 82. As- 
suming a population of red giants and applying the 
diagnostics of Origlia et al. (1997), the EWs imply 
[Fe/H] « —0.5 to —0.2 dex for a carbon depletion 
of [C/Fe] = 0.0 to —0.5 dex. Larger carbon deple- 
tions are ruled out since for [C/Fe] < —0.5 dex, the 
OH bands at 1.6265 /*m become comparably deep 
or deeper than the 12 CO (6,3) bandhcad (Origlia 
et al. 1997), which is not the case in the 3D spec- 
tra. Hence, the 3D data are definitely inconsis- 
tent with a dominant old metal-rich population in 
the central 35 pc at the nucleus. Metallicity esti- 
mates for supergiants using the diagnostics from 
Oliva & Origlia (1998) depend more sensitively on 
T e g and are less well constrained. For the derived 
T off « 3800 K, [Fe/H] « -0.5 to 0.0 dex depend- 
ing on [C/Fe]. Lower temperatures would reduce 
[Fe/H] by up to 0.3 dex while higher temperatures 
up to 4300 K would increase it by up to 0.4 dex. 
The analysis of the EWs in section 5.1.2 leading to 
a consistent interpretation together with the near- 
solar abundances of the H II regions support that 
the supergiants in the central 35 pc of M 82 have 
roughly solar metallicity. 

Similar metallicities are inferred over the en- 
tire regions mapped with 3D, assuming a domi- 
nant population of supergiants and accounting for 
the variations in T c g. The spatial variations in 
the CO bandhead EWs are not likely due to vari- 
ations in [Fc/H] of the stars. Indeed, the observed 
ranges for Wi.62 and W2.29 would imply variations 
in the metallicity by factors of at least 4 — 5, not 
plausible on scales of ~ 10 — 100 pc and over the 
10 — 50 Myr lifetimes of red supergiants. 

5.2. Additional Continuum Emission Sources 

The possible sources responsible for the dilu- 
tion of the stellar absorption features in M 82 in- 
clude young OB stars, nebular free-free and free- 
bound processes, and dust heated at 600 — 1000 K 
by the OB stars ("hot dust"). We estimated 
the broad-band emission from OB stars using the 
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number of representative 08.5 V stars given in ta- 
ble 9 and the photometric properties tabulated by 
Vacca et al. (1996) and Koornnecf (1983). We 
computed the contribution of the nebular con- 
tinuum emission from the dereddened Br7 fluxes 
using the relationships given by Satyapal et al. 
(1995). These are for case B recombination with 
n c = 100 cm" 2 and T c = 10 4 K. They are little 
affected by the electron density, but depend more 
sensitively on the electron temperature. However, 
for T c « 5000 K, the nebular flux densities inferred 
from Br7 would be lower (e.g. Joy & Lester 1988). 

For selected regions (table 10) and across the 
entire 3D field of view, OB stars and nebular pro- 
cesses make a negligible contribution to the near- 
infrared continuum emission (see also Satyapal et 
al. 1995), leaving the hot dust as most important 
source of dilution. A crude estimate of the di- 
lution for the 3D field of view can be obtained 
independently from the SWS data. We assumed 
grey body emission for the hot dust and adopted 
?hd = 800 K and n = 1.5, consistent with previ- 
ous work (e.g. Smith et al. 1990; Larkin et al. 
1994), and also with the negligible dilution near 
1.6 /itm. With extinction-corrected continuum flux 
densities of ~ 9 x 10~ 13 Wm _2 /im _1 near 4^tm 
and « 3 x 10~ 12 Wm _2 ^m _1 near 2.3 ^m, the 
predicted contribution from hot dust at the lat- 
ter wavelength is 10% — 15%, consistent with the 
dilution inferred from the EWs alone. 

5.3. Extinction towards the Evolved Stars 

We finally constrained the extinction towards 
the evolved stellar population, important for de- 
riving its intrinsic luminosity since it has a very 
different spatial distribution than the ionized gas 
(see figure 3) and may suffer from different levels 
of obscuration (see also e.g. McLeod et al. 1993). 

The extinction for selected regions was derived 
from minimum x 2 -htting to the 3D spectra as fol- 
lows. For each region, we combined the JT-band 
spectrum for the appropriate stellar spectral type 
with a grey-body emission curve for the hot dust 
with Thd = 800 K and n = 1.5, in the propor- 
tions given by -D2.3, and adjusted the extinction 
(Ay) for the best fit to the observed spectrum. 
We considered a uniform foreground screen and 
a mixed model, and adopted the extinction law 
from Draine (1989). The template stellar spectra 
were taken from the atlases of Forster Schrciber 



(2000) and Kleinmann & Hall (1986), convolved 
to the spectral resolution of the M82 data when 
appropriate. Due to the rather poor sampling for 
K supcrgiants, a K5 I template was used for all 
regions except Bl, for which the K0 I and K5 I 
spectra available were averaged to produce a tem- 
plate K2 I spectrum. We excluded the _ff-band 
data from this analysis because available libraries 
had too limited a wavelength coverage. 

The results are not significantly affected by the 
choice of template star within a few spectral sub- 
classes, by Thd in the range 600 — 1000 K, by n 
between 1 and 2, and by the power-law index for 
the extinction law within 0.1 — 0.2 dcx. The as- 
sumption that the same extinction applies to the 
evolved stars and to the hot dust is of little con- 
sequences since the stars dominate the continuum 
emission. We have also constrained the extinction 
from the H — K colour excess (after correction 
for dilution), using the stellar data compiled by 
Koornneef (1983). The values obtained are in ex- 
cellent agreement with those from the spectral fits. 

The extinction towards the evolved stars is 
lower than towards the ionized gas, except for the 
central 35 pc. The extinction corrections differ 
by up to ss 35% for purely foreground obscura- 
tion and up to ~ 60% for the mixed model, pre- 
sumably reflecting the different distributions for 
the corresponding sources relative to the obscur- 
ing dust. The data do not allow the distinction be- 
tween different model geometries owing to the rel- 
atively small wavelength coverage. However, the 
extinction corrections are similar to within 35% 
or less, except for the central 35 pc (« 50%). We 
will adopt the results for the mixed model because 
such a geometry seems more plausible for popula- 
tions of clusters of supcrgiants, which are likely 
more or less uniformly mixed with interstellar gas 
and dust clouds, as for the H II regions. 

We derived the extinction across the 3D field 
of view from the dilution-corrected H — K colour 
map as described above, for rebinned 1" x 1" pix- 
els and assuming an intrinsic population of super- 
giants of appropriate type at each location. We 
considered only foreground obscuration along any 
line of sight. Comparing the flux density inte- 
grated over the rebinned, dilution- and extinction- 
corrected /i'-band map to the observed integrated 
flux density, the global correction factor at 2.2 /im 
for this non-uniform foreground screen model is 
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very close to that obtained from the fits to the 
3D field of view spectrum; it corresponds to an 
effective Af IX = 18 mag (or A^ FS = 8 mag). 

Figure 11 shows the results of the popula- 
tion/spectral synthesis for selected regions and 
also compares the 3D iJ-band spectra corrected 
for extinction to the template K4 I spectrum from 
Dallicr et al. (1996) — their only K-type super- 
giant for solar metallicity. This figure illustrates 
well the quality of the fits, in particular for the 
numerous stellar absorption features. The num- 
ber of representive stars required to produce the 
dilution- and extinction-corrected Lk, computed 
using the stellar data compiled by Schmidt-Kalcr 
(1982) and Koornneef (1983), are reported in ta- 
ble 10. The table also gives the properties for the 
starburst core which are derived in appendix B to- 
gether with additional constraints (bolometric lu- 
minosity, mass, and rate of supernova explosions) 
used for the starburst modeling of paper 2. 

6. SUMMARY AND DISCUSSION 

We have obtained near-infrared imaging spec- 
troscopy with the 3D instrument and mid-infrared 
spectroscopy with the ISO-SWS of the starburst 
regions of M 82. We have used these data, together 
with results from the literature, to determine the 
physical conditions of the ISM and the compo- 
sition of the stellar population on spatial scales 
ranging from a few tens of parsecs to 500 pc. 

The central regions of M 82 show complex and 
important structure on scales at least as small as 
~ 25 pc, with the tracers of H II regions and those 
of cool evolved stars having in addition very differ- 
ent spatial distributions. However, the structural 
properties within the most recent star-forming re- 
gions and the excitation state of the nebular gas 
are remarkably similar throughout the entire star- 
burst core, on all spatial scales. The picture which 
emerges is that of a large number of closely-packed 
ionizing stellar clusters and small gas clouds, with 
a highly homogcneized distribution within larger- 
scale concentrations. The near constancy in local 
structural properties (as reflected notably in the 
ionization parameter) and in dominant OB star 
population suggests a similar star formation effi- 
ciency and evolutionary stage for the most recent 
star-forming regions over the entire central 500 pc. 

The variations in the near-infrared continuum 



properties are partly attributable to spatially non- 
uniform extinction and contribution from hot 
dust emission, but also to variations in the in- 
trinsic composition of the evolved stellar popula- 
tion. Together with the M*/Lk ratio, the stel- 
lar absorption features are consistent with solar- 
metallicity red supergiants as main sources of 
near-infrared continuum emission throughout the 
starburst core, down to a few tens of parsecs at 
the nucleus. Although the derived spectral types 
cover a fairly large range, they correspond in fact 
to a rather narrow range in evolutionary stages. 

The OB stars and red supergiants, which dom- 
inate the luminosity of M82, thus trace the star 
formation history within the starburst core up to 
about 50 Myr ago. The complex starburst history 
in M82 is already obvious from the relative dis- 
tributions of the ionized gas and of the red super- 
giants, suggesting a time sequence in the triggering 
of the bursts at different locations. Radial evolu- 
tion has been proposed in the past but some au- 
thors favoured inside-out propagation while others 
favoured the opposite scenario (e.g. McLcod et 
al. 1993; Shcn & Lo 1995; Satyapal et al. 1997; 
de Grijs et al. 2000). We will address this issue in 
paper 2 by applying starburst models to the data 
presented here in order to constrain quantitatively 
the detailed spatial and chronological evolution of 
starburst activity in M82. 
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A. DERIVATION OF THE EFFECTIVE IONIZATION PARAMETER IN M 82 



In this appendix, we present our determination of the effective ionization parameter for the photoionized 
nebulae in M82. Since it depends on the flux of Lyman continuum photons impinging on the gas (see 
equation [5], section 4.3), its derivation requires a detailed knowledge of the properties and distribution 
of the ionizing stars and gas clouds. We selected three representative regions for this analysis: the entire 
starburst core, and regions Bl and B2. For clarity, most numerical results are reported directly in table 11. 

A.l. Properties and Distribution of the Gas Clouds 

Various observations of the molecular, neutral atomic, and ionized gas in the starburst regions of M 82 
reveal important structure on scales at least as small as w 20 — 30 pc (e.g. figure 3; Larkin et al. 1994; Shcn 
& Lo 1995; Achtermann & Lacy 1995; Satyapal et al. 1995). Models of the ISM imply even more extreme 
properties, suggesting that the gas clouds possess small and dense molecular cores with thin neutral atomic 
surfaces, and are embedded in large ionized envelopes in a highly pressurized ISM (Olofsson & Rydbeck 
1984; Lugtcn et al. 1986; Duffy et al. 1987; Wolfire, Tielens, & Hollenbach 1990; Lord et al. 1996; Stutzki 
et al. 1997). We adopted the neutral cloud properties derived by Lord et al. (1996). For the starburst core, 
we took their average cloud radius r c \ and mass M c \ inferred from the global properties of M 82. For Bl and 
B2, we adopted their results for the southwestern infrared emission lobe encompassing these regions. 

Combining M c \, the mass of molecular gas Mh 2 , and the volume V in which the clouds are distributed 
yields their space number density n c \ and mean separation <i c i_ c i. Since V is not well constrained, we 
considered two limiting cases for each region. For the starburst core, we adopted a sphere of diameter 500 pc 
(i.e. 30"), and an edge-on disk of radius and thickness of 200 pc more consistent with the global distribution 
of various gas components (e.g. Larkin et al. 1994; Shen & Lo 1995; Achtermann & Lacy 1995). For Bl 
and B2, we assumed spheres of radius 19.5 pc, and columns with radii of 19.5 pc and lengths equal to the 
intersection along the line of sight of the starburst disk at the corresponding projected locations (230 pc and 
365 pc, respectively). The cross-section radii were chosen to cover an area equivalent to the 2.25" x 2.25" 
square aperture used to extract the spectra from the 3D data cubes. 

We estimated Mh 2 at Bl and B2 from the CO J = 1 — > map of Shen & Lo (1995), which has a spatial 
resolution comparable to the 3D data (2.5"). We converted the CO intensities into H2 column densities 
and masses using Nn 2 /Ico 1^0 = 7 x 10 19 cm _2 K _1 (kms -1 ) 1 derived by Wild et al. (1992) for regions 
encompassing Bl and B2. For the entire starburst core, these authors estimated Mh 2 = 1.8 x 10 8 M , the 
value adopted by Lord et al. (1996). We emphasize that Wild et al. derived the A/h 2 /-^CO i-»o conversion 
factor from detailed radiative transfer calculations applied to observations of 12 CO (up to J = 6 — > 5), 13 CO, 
and 18 CO lines at various positions along the galactic plane of M 82, so that it should account properly for the 
molecular gas mass at each location. We also derived the mass of ionized gas Mjj+ from the Ql yc determined 
in section 4.1, assuming case B recombination with appropriate electron density and temperature for M82 
(see e.g. Osterbrock 1989 for details). The radius n of the outer ionized surface of the clouds is set by 

where N c i is the total number of clouds and $y + is the volume filling factor of the ionized gas. The n does 
not depend on the volume considered because does not. 

A. 2. Properties and Distribution of the Ionizing Stellar Clusters 

The OB stars in M82 are likely to reside mainly in clusters. Indeed, high-resolution HST imaging of 
the central regions of M82 reveals the presence of over a hundred compact, luminous clusters with ages 
estimated between ~ 10 Myr and 10 Gyr (O'Connell et al. 1995; Gallagher & Smith 1999; de Grijs et al. 
2000). Due to the high extinction in M82, these are probably mostly foreground, but some may reside in 
the inner regions if they lie in directions of lower extinction. Near-infrared observations provide evidence for 
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the presence of such systems deep in the obscured nuclear regions of M82: several distinct compact if-band 
continuum sources are seen, with surface brightnesses, sizes, and CO bandhead strengths consistent with 
young clusters of red supergiants (Satyapal et al. 1997; see also section 5). Such super star clusters are 
observed in a growing number of starburst systems and may constitute an important mode of star formation 
in starbursts (e.g. Holtzmann et al. 1992; O'Connell, Gallagher, & Hunter 1994; Whitmore & Schweizer 
1995; Maoz et al. 1996; Ho & Filippenko 1996; Tacconi-Garman, Sternberg, & Eckart 1996). 

The space number density n+ and mean separation of the ionizing clusters can be computed from 
the intrinsic Lyman continuum photon emission rates. We assumed that the clusters follow the luminosity 
function (LF) 

acfe-WW^ (A2) 

with (3 = 0.19 in the range Q£ yc = 10 45 - 10 49 5 s" 1 and (3 = 1 in the range Ql yc = 10 49 5 - 10 53 s" 1 . 
This LF is derived by Thornlcy et al. (2000) based on the optical and Lyman continuum LF's observed 
for super star clusters and H II regions in a variety of local star-forming galaxies, with extension to low 
luminosities from Monte-Carlo simulations. The upper limit in Q^ yc is consistent with the properties of the 
optically-selected clusters in M82. O'Connell et al. (1995) measured intrinsic absolute U-band magnitudes 
in the range —9.6 mag to —14.5 mag. For an unevolved stellar population with a Salpeter (1955) initial 
mass function, these values correspond to Q^ yc = 10 50 5 — 10 52 ' 7 s _1 . Accounting for possible evolutionary 
effects, the U-band magnitudes are still consistent with high Q^yc U P to ~ 10 52 s _1 . The lower limit of the 
LF corresponds to the smallest associations susceptible of ionizing an H II region, i.e. containing at least 
one early-B star. 

The results, summarized in table 11, suggest that the starburst regions of M82 can be represented by a 
collection of closely-packed ionizing clusters and small gas clouds, separated on average by a few parsecs. The 
comparable mean separations inferred for the clouds and for the clusters suggest that they have a well-mixed 
distribution throughout the regions of interest, as illustrated in the bottom of figure 7. Uncertainties in the 
LF or in the cloud properties are not likely to alter this picture. Substantial differences in n± and require 
cither a very flat LF ((3 < 0.5 at high luminosities), or a very high lower cutoff (near the inflection point at 
Qlyc = 10 49 ' 5 s _1 ). The former case would not be consistent with the observed range (3 = 0.5 — 1.0 (see 
Thornley et al. 2000 and references therein), while the latter would imply the unlikely situation that only 
massive clusters containing at least one 50 M star can form, excluding smaller OB associations. Similarly, 
the cloud properties would need to be very different than those adopted. For example, increasing the M c \ 
and r c i for Bl and B2 to the values for the starburst core and vice- versa implies a variation by 36% only for 

A. 3. Effective Ionization Parameter 

Adopting therefore a random distribution model, we derived U e s by considering the Lyman continuum 
photon flux incident on the surface area afforded by the outer ionized edge of the clouds. By analogy with 
equation (5), 

^Lyc 

4irrfN cl n c c 

The resulting log[/ c ff are —2.3 to —2.4 dcx. For comparison, we also computed values of \ogU applying 
directly the central cluster model with the radius of the spherical volumes considered above. The U e g values 
are lower by about an order of magnitude, due to the important increase in surface area of the gas exposed 
to the Lyman continuum radiation field for the more realistic randomized distribution. 

Variations in the input parameters over plausible ranges do not affect importantly the derived logJ7 e ff. 
For instance, varying n e in the range 10 2 — 10 3 cm~ 3 changes the thickness of the ionized layers but the effect 
on log U e H due to the variation in exposed surface is partly compensated by that of n c itself, resulting in 
differences smaller than ±0.2 dex. A higher T e = 10 4 K decreases log U c g by about 0.15 dex. Interverting the 
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neutral cloud properties between the small-scale regions Bl and B2, and the starburst core implies \ogU c g 
higher and lower by about 0.1 dex, respectively. The logt/ e ff is fairly well constrained in the range —2 to 
—2.6 dex; we have adopted a representative —2.3 dex. 

Interestingly, very similar conditions are derived for individual regions on scales of a few tens of parsecs 
as well as for the 500-pc size starburst core. Bl and B2 contribute however only about 5% of the total <3L yc 
determined for the entire starburst core. This indicates little variation in the degree of ionization of the 
photoionized nebulae throughout M82. In order to better understand this result physically, equation (A3) 
can be expressed in terms of more fundamental properties. Neglecting the thin neutral atomic layer, N c \ = 
Mh 2 /M c i and M c \ = 7rr^ 1 A/'cimH, where Af c i is the H column density of a molecular cloud core and tor is the 
mass of a H atom. The mean free path for Lyman continuum photons is determined by the neutral surface of 
the clouds, AL yc = (n c i7rr^j) _1 . We also introduce A, = (n c \irrf ) _1 which can be interpreted as a geometrical 
mean free path characterizing the path required by a Lyman continuum photon to reach the outer ionized 
surface of a cloud. Combining these relationships, equation (A3) can be re-written as 

The effective ionization parameter is thus related to Qhyc/Mn 2 which provides a measure of the star formation 
efficiency. Consequently, for similar gas cloud properties, the comparable U e s at different locations and on 
various spatial scales results from a near constancy in the average star formation efficiency throughout the 
starburst regions of M82. 



B. ADDITIONAL OBSERVATIONAL CONSTRAINTS FOR M 82 

In the following, we derive further constraints that will be used in paper 2 for the application of starburst 
models to M 82. These include the near-infrared continuum properties for the starburst core, as well as the 
bolomctric luminosity, the stellar mass, and the rate of supernova explosions for various regions. 

B.l. Near-infrared Properties of the Starburst Core 

Satyapal et al. (1997) measured the depth of the CO bandheads longwards of 2.3 ^m in a 24"-diameter 
aperture on the nucleus of M82. From a calibration derived using the stellar atlases of Kleinmann & 
Hall (1986) and Forster Schreiber (2000), their spectroscopic index of 0.18 mag (which is independent of 
extinction) corresponds to W2.29 — 14.5 A. This is essentially the global value measured for the 3D field of 
view accounting for 10% dilution (14.4 A). Dilution can be neglected for the starburst core, since about half 
of the infrared emission by dust and of the nebular emission from H II regions originates outside of the 3D 
field of view while this fraction is 80% for the if-band emission (see table 10 and section 4). Because of 
the similar intrinsic W2.29 for the 3D field of view and the central 24" as well as the relative constancy of 
M*/Lk within the central 500 pc, we adopted the W1.62 for the 3D field of view as representative for the 
starburst core. 

While the observed if-band magnitudes determined by various authors for the central w 500 pc of M 82 
agree well with each other, there are substantial discrepancies in the reported absolute intrinsic magnitudes 
(M&). The results vary between -23.3 mag and -22.0 mag (Rieke et al. 1980; Telesco et al. 1991; McLeod et 
al. 1993; Satyapal et al. 1997). These differences are mainly attributable to the different extinction values and 
model geometries assumed. We re-examined this issue using the results from our 3D and SWS observations. 
Based on the morphology of the if-band and nebular line emission and on the global distributions of the 
corresponding sources, we argue as Rieke et al. (1980) that half of the extinction towards the bulk of ionized 
gas probably applies to the bulk of evolved stars. Taking the observed -ftf-band magnitude measured by 
Telesco et al. (1991, 5.43 mag), an Ay lx — 26 mag implies = —23.2 mag. Similar values are obtained 
with the mixed extinction model for the 3D field of view (—22.9 mag) and the effective extinction derived 
from the H — K map (—23.0 mag). Our results are consistent with the lower range of values reported 
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in the literature. The 3D data trace projected regions close to the nucleus which are on average the most 
obscured, leading to a possible overestimate of the global extinction for the starburst core. On the other hand, 
large optical depths are inferred, which may indicate that the near-infrared light does not probe the stars 
throughout the entire galaxy along the line of sight. Therefore, we will adopt an average = —23.0 mag, 
implying an effective global Ay lx — 21 mag (or A V FS = 9 mag) towards the evolved stars in M82. 

B.2. Bolometric Luminosity 

Two main sources dominate the bolometric luminosity (Lboi) in M82: the hot massive stars and the cool 
evolved stars. In dusty starbursts such as M82, the infrared luminosity (Lir) provides a good approximation 
to the Lboi from hot stars. For the entire starburst core, we adopted Lir = 3x 10 10 L Q (Telesco & Harper 
1980). We estimated Lir for individual regions using measurements of the mid-infrared emission which has 
a similar morphology as the infrared emission out to at least 100 /im, where the global energy distribution 
of M82 peaks (Telesco et al. 1991). Telesco et al. (1993) obtained Lir = 18 Ljy for the starburst core of 
M 82, where Ln is the luminosity in the TV-band centered at 10.8 /im, also valid within a factor of two for a 
sample of 11 starburst galaxies. We actually used the 12.4 ^m map of Telesco & Gezari (1992) 7 which has 
a spatial resolution comparable to the 3D data, scaling the fluxes according to /jv = 0.4/i2.4 Mm derived by 
these authors. Combining the above relationships, 

^ = 3.25 x 10 7 ( ( f -^-) . (Bl) 
L© \MpcJ \ Jy J 

Equation (Bl) includes average extinction effects as well as contributions in the 12.4 jum and iV-bandpasses 
from PAH emission features, emission lines, and the silicate absorption feature. From our SWS and 3D data 
together with the fine-structure line maps of Achtermann & Lacy (1995) and the 3.3 /zm PAH map of Satyapal 
et al. (1995), we estimate that spatial variations in the spectral features and in extinction introduce errors 
of about 35% in applying equation (Bl) to individual regions. Following McLeod et al. (1993), we included 
an additional 30% of the detected Lir in Lboi to account for light escaping in directions perpendicular to 
the galactic plane. 

Individual giants and supergiants with temperatures between 3500 K and 6000 K have L^ \/Lk ~ 10 — 30, 
similar to the ratio for mixed populations of cool evolved stars (e.g. McLeod et al. 1993). We have thus 
estimated the Lboi from evolved stars using the intrinsic stellar Lk and assuming a representative ratio of 
20, with 50% uncertainty Table 12 summarizes the results for selected regions. 

B.3. Mass 

We estimated the dynamical mass M& yn from position-velocity maps available in the literature. These 
include observations of the CO J = 1 — > millimeter line (Shen & Lo 1995), of the [Ne II] 12.8 ^m line 
(Achtermann & Lacy 1995), and of the [S III] A9069 A line (McKeith et al. 1993). We assumed a uniform 
mass distribution and dynamical equilibrium, so that within projected radius r 

M dy „« r) = 3.49xl0»(^)(^) 2 , <B2) 

where v rot is the rotational velocity. Since various observations of the molecular and ionized gas reveal 
important concentrations which are interpreted as circumnuclcar rotating rings or spiral arms (e.g. Larkin 
et al. 1994; Shen & Lo 1995; Achtermann & Lacy 1995; Seaquist et al. 1998; Neiningcr et al. 1998), we 
determined v rot using the peak velocity at the locations of these concentrations, corrected for an inclination 
of 80°, and for beam and velocity smearing. In addition, we used the dynamical mass obtained by Gaffney et 
al. (1993) in the central 1" of M82, from stellar velocity dispersion measurements using the CO bandhead 

kindly made available to us in digital form by the authors 
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at 2.29 /im. The resulting Md yn versus r curve is plotted in figure 12 along with the mass model proposed 
by Gotz et al. (1990). Our results imply larger masses for r < 10", by up to a factor of « 4 at r = 0.5". The 
more recent data used here have higher resolution than those used by Gotz et al. (1990), and thus probe 
much better the central regions of M82. 

To obtain the stellar mass M*, we subtracted the gaseous mass from Md yn - The mass of ionized and 
molecular hydrogen (Mh+ and Mh 2 ) for the starburst core is estimated in appendix A. Following the same 
procedure as described there for Bl and B2, we computed M H + and Mjj 2 for the central 35 pc from the 
intrinsic Qhyc, and from the CO data of Shen & Lo (1995) with the conversion factor between the CO 
intensity and H 2 column density from Wild et al. (1992) at the position of the nucleus. Table 13 gives the 
various mass estimates. The resulting M* represent lower limits since the Mh 2 adopted for the starburst 
core was determined by Wild et al. (1992) in a region larger than 30" in diameter, and the gas observed 
towards the central 35 pc of M 82 is probably mainly located in a circumnuclear ring at larger radius, as 
mentioned above. 

B.4. Rate of Supernova Explosions 

Estimates of the global rate of supernova explosions (i*jn) m M82 from the properties of the compact, 
non-thermal synchrotron emission sources detected at centimeter wavelengths (size and radio luminosity 
distributions, luminosity variations) vary in the range 0.02 — 0.1 yr _1 (e.g. Kronberg & Sramek 1985; van 
Buren & Greenhouse 1994; Huang et al. 1994; Muxlow et al. 1994; Allen & Kronberg 1998). For selected 
regions in M82, we applied the relationship between j^sn and [Fe II] 1.644 /im line flux derived by Vanzi 
& Rieke (1997). Our fluxes corrected for the extinction inferred from the Brackett lines imply rates of 
- 3 x 10~ 3 yr" 1 for the central 35 pc of M82 and Bl, ~ 6 x 10~ 3 yr" 1 for B2, and ~ 0.01 yr" 1 for the 
3D field of view. We emphasize however that the uncertainties in i'sn's inferred from [Fe II] fluxes are 
very large. The calibration of Vanzi & Rieke (1997) is based on the integrated [Fe II] luminosity and radio 
vsn, and depends on uncertain assumptions about the supernova remnants (SNRs) lifetimes. The [Fe II] 
emission may trace a different population of SNRs than radio observations, and its interpretation may be 
further complicated by other excitation sources such as an outflowing starburst wind (e.g. Greenhouse et 
al. 1997). A calibration based on the average intrinsic line fluxes of four compact [Fe II] sources candidate 
SNRs measured by Greenhouse et al. (1997) would imply ^sn's about 50 times smaller! 
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Fig. 1 — Regions observed in M82. The 3D field 
of view and ISO-SWS apertures are shown on a K- 
band map (Forster et al. 1994) with superimposed 
contours of the CO J = 1 — > emission (Shen & 
Lo 1995). The cross and the triangle indicate the 
positions of the nucleus and of the western mid- 
infrared emission peak, respectively. 

Fig. 2. — 3D H- and if-band spectra of selected 
regions in M82. The top three panels show the 
spectra of the nucleus and of regions Bl and B2 
taken in 2.25" x 2.25" apertures (see table 2) while 
the bottom panel shows the spectrum of the en- 
tire 3D field of view. The effective resolution is 
R <~ 1015 in the iJ-band and R ~ 830 in the K- 
band. The vertical axis is a linear flux density 
scale. The spectra are normalized to unity in the 
interval 2.2875 - 2.2910 /im. The absolute flux 
densities in Wm~ 2 /im _1 can be recovered using 
the multiplicative factors 6.52 x 10~ 14 for the nu- 
cleus, 2.96 x 10~ 14 for B2, 1.78 x 10~ 14 for Bl, and 
7.24 x 10" 13 for the 3D field of view. The posi- 
tions of various lines are indicated on the spectra 
of the nucleus (stellar absorption features) and of 
Bl (emission lines). 

Fig. 3. — Selected near-infrared maps of M 82 ob- 
tained with 3D. The greyscale levels are indicated 
next to each panel, in the units given in the fol- 
lowing, (a): K-hand emission map, with con- 
tours of the -ff-band emission from 0.5 to 2.15 
in steps of 0.15; greyscales and contours arc in 
units of 10~ 14 Wm~ 2 /iin^ 1 arcsec~ 2 . (b): Br7 
emission line map, with contours of the He I 
2.058 fim line emission from 0.5 to 1.5 in steps 
of 0.1; greyscales and contours arc in units of 
10~ 17 Wm~ 2 arcscc~ 2 . (c): Map of the equiva- 
lent width of the 12 CO (2,0) bandhead at 2.29 fim 
(W2.29) in units of A; contours of the equivalent 
width of the 12 CO (6,3) bandhead at 1.62 /im 
(W1.62) from 4.8 A to 6.6 A, in steps of 0.3 A. 
(d): log(L K /L Lyc ) map. (e): He I 2.058/Br 7 line 
ratio map. (f): log(W/i. 62/^2. 29) map. The axis 
coordinates are relative offsets from the nucleus, 
indicated by the cross. From left to right, the 
boxes in each panel indicate the apertures used to 
extract the spectra of the central 35 pc of M82, 
B2, and Bl plotted in figure 2 (labeled "N," "B2," 
and "Bl" in panel a). 



Fig. 4. — ISO-SWS mid-infrared spectrum of M 82 
(full scan AOT SWS01). The spectral resolution 
varies from ~ 1000 at short wavelengths to ~ 500 
at long wavelengths. The "jumps" in the contin- 
uum level at 12.0 fim, 27.8 /im, and 29.5 /mi arc 
caused by the increase in aperture size and the 
fluctuations in the continuum, especially in the 
4 — 5 /im region, are mainly due to noise. 

Fig. 5. — High S/N ratio spectra of individual 
mid-infrared lines in M 82, obtained with the ISO- 
SWS (line scans AOT SWS02). The spectral res- 
olution ranges from ~ 2000 to ~ 1000 from short 
to long wavelengths. 

Fig. 6. — Results of the determination of the 
global extinction towards the ionized gas in M82. 
The plots show the intrinsic Lyman continuum 
photon emission rates Ql yc for the SWS 14" x 20" 
aperture, derived from the H recombination line 
fluxes corrected for the best-fit extinction for dif- 
ferent geometries as indicated in each plot. Pan- 
els to the left show the results using the Draine 
(1989) extinction law throughout the entire in- 
frared range while those to the right show the 
results assuming the extinction law towards the 
Galactic Center between 3 /xm and 10 /im given 
by Lutz (1999). The horizontal lines in each plot 
indicate the average <5L yc . 

Fig. 7. — Geometries considered in the derivation 
of the effective ionization parameter in M 82. Top: 
a gas shell of radius R surrounding a central stel- 
lar cluster, or association of clusters. Bottom: a 
well-mixed distribution of stellar clusters and gas 
clouds. 

Fig. 8. — T e g-sensitive ratios of mid-infrared 
fine-structure lines (top three panels) and of near- 
infrared recombination lines (bottom two panels). 
The various curves represent the theoretical pre- 
dictions for different sets of parameters. The solid 
lines show the results for the nebular parameters 
derived for M82: n-n ~ n c — 300 cm~ 3 , logU = 
—2.3 dex, R — 25 pc, solar gas-phase abundances, 
and no interstellar dust grains mixed with the ion- 
ized gas. The effects of changing the ISM compo- 
sition to a gas and dust mixture as in the Orion 
nebula, the gas density to rin = 10 3 cm~ 3 , or log U 
between —2 dex and —2.5 dex are illustrated as 
well (see labels in each plot). The horizontal bars 
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show the extinction-corrected ratios obtained with 
the SWS and with 3D, with vertical width indi- 
cating the measurement uncertainties. Different 
shading or filling patterns in the two bottom pan- 
els show the ratios for selected regions: the central 
35 pc at the nucleus (empty bar and horizontal line 
marking the upper limit, labeled "N"), Bl (dark- 
shaded bar), B2 (cross-hatched bar), and the 3D 
field of view (light-shaded bar) . 

Fig. 9. — Spectroscopic indices sensitive to the ef- 
fective temperature and luminosity class of cool 
stars. The data for the central 35 pc of M82 
(labeled "N"), Bl, and B2 are indicated by the 
horizontal bars, with vertical width corresponding 
to the measurement uncertainties. The integrated 
equivalent widths for the 3D field of view are simi- 
lar to those of B2. The stellar data are taken from 
the compilation by Forster Schrcibcr (2000); open 
circles, filled circles, and crosses represent super- 
giants, giants, and dwarfs, respectively. 

Fig. 10. — Diagnostic diagrams for the estimation 
of the amount of dilution near 1.6 /im and 2.3 /xm. 
The shaded areas indicate the loci of giants (dark 
shade) and supergiants (light shade), based on the 
stellar data compiled by Forster Schreiber (2000). 
The arrows indicate the effects of dilution by fea- 
tureless continuum sources. Stellar effective tem- 
peratures corresponding to various intrinsic Wi.62 
are labeled on the right-hand side diagrams. Top: 
data for the central 35 pc of M82 (labeled "N"), 
Bl, B2, and the 3D field of view (labeled "3D"). 
Bottom: data for individual 1" x 1" pixels from 
the rebinned 3D maps, with typical uncertainties 
shown in the upper left corner of the diagrams. 

Fig. 11. — Results of spectral synthesis for selected 
regions in M82. The results of the fits are shown 
for the central 35 pc of M82 (top panels), regions 
Bl and B2 (middle panels), and the 3D field of 
view (bottom panels). The black lines arc the 
3D spectra corrected for the best-fit mixed extinc- 
tion (table 10). The grey lines are the combina- 
tion of the template stellar spectrum and the hot 
dust ("HD") emission appropriate for each region. 
The _ff-band template spectrum is the K4 I star 
BS5645 from Dallier et al. (1996). The if-band 
template spectra, from Kleinmann & Hall (1986), 
are the K5 I star BS 8726, and a K2 I obtained by 
averaging the spectra of BS 8726 and of the K0 I 



star RW Cep. 

Fig. 12. — Enclosed dynamical mass versus pro- 
jected radius in M 82 (derived in appendix B). The 
diagonal line represents the mass model proposed 
by Gotz et al. (1990). The horizontal line seg- 
ments at projected radii of 1.3" and 15" indicate 
the contribution from the ionized and molecular 
hydrogen gas for the central 35 pc and the entire 
starburst core. 
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Table 1 

Log of the 3D observations of M 82 



Date 


Telescope a 


Field b 


Band 


tint ° 

(s) 


Seeing 
(arcsec) 


Atmospheric calibrator 


1995 Jan 13 


CA 


1 


K 


600 


1.5 


HD 82189 (F5 V) 






2 


K 


600 


1.5 


HD 82189 (F5 V) 


1995 Jan 14 


CA 


3 


K 


600 


1 


PPM 17105 (GO V) 






3 


H 


480 


1 


HD 87141 (F5 V) 


1995 Jan 16 


CA 


4 


K 


600 


1.3 


PPM 17105 (GO V) 


1995 Jan 21 


CA 


1 


H 


480 


1 


HD 87141 (F5 V) 






4 


H 


480 


1 


HD 87141 (F5 V) 


1996 Jan 06 


WHT 


2 


H 


900 


1 


HD 26356 (B5 V) 



a CA : 3.5 m telescope at Calar Alto, Spain. WHT : 4.2 m William-Herschel- Telescope 
on the Canary Islands, Spain. 

b The right ascension and declination offsets of the fields with respect to the nu- 
clear position (defined by the if -band emission peak at 0:1950: 09 h 51 m 43?53, 61950- 
+69°55'00'.'7; Dietz et al. 1986) are the following: field 1, (+2", +1.5"): field 2, 
(+2", -2"); field 3, (-3",-l"); field 4, (-8", -3"). 

c Total on-source integration time per detector pixel. 
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Table 2 

Observed near-infrared properties of selected regions in M 82 



Band or feature Nucleus Bl B2 3D field 

(0",0") a (-f0",-4.25") a (-5.25", -2") a 

Broad-band flux densities b ( Jy) 



ii-band 


0.076 


0.021 


0.032 




1.01 


if -band 


0.104 


0.030 


0.048 




1.20 


Emission line fluxes 


(10" iY Wm" 2 ) 








[Fe II] a 4 F 9/2 - a 4 D 5/2 (1.5335 /xm) 


< 0.11 


0.75 ±0.09 


0.56 ±0. 


.04 


9.2 ±2.9 


H I Brl3 (n = 4 - 13, 1.6109 /xm) 


0.48 ±0.30 


0.90 ±0.08 


0.71 ±0. 


.16 


13.2 ±4.5 


H I Brl2 (n = 4- 12, 1.6407 (im) 


0.33 ±0.08 


0.81 ±0.19 


0.75 ±0. 


45 


12.9 ±6.5 


[Fe II] a 4 F 9/2 - a 4 D 7/2 (1.6435 /xm) 


4.12 ±0.07 


3.53 ±0.44 


4.99 ±0. 


44 


137 ±12 


[Fe II] a 4 i^ 7/2 - a 4 D 5/2 (1.6769 /xm) 


0.52 ±0.37 


0.22 ±0.07 


0.34 ±0. 


.13 


10.5 ±3.6 


H I Brll (n = 4- 11, 1.6807 /xm) 


0.42 ±0.25 


1.44 ±0.28 


0.90 ±0. 


19 


14.3 ±4.1 


He I 3 3 P - 4 3 L> (1.7008 /xm) 


< 0.16 


0.50 ±0.12 


0.38 ±0. 


09 


6.3 ±2.8 


H I BrlO (n = 4 - 10, 1.7362 /xm) 


1.22 ±0.24 


2.31 ±0.18 


1.76 ±0. 


17 


32.2 ±5.2 


He I 2 1 S-2 1 P (2.0581 /xm) 


2.53 ±0.28 


5.19 ±0.08 


4.34 ±0. 


14 


74.9 ±3.3 


H 2 1-0 S(l) (2.1213 /xm) 


1.01 ±0.21 


0.66 ±0.07 


0.98 ±0. 


11 


21.3 ±2.4 


H I Br 7 (n = 4 - 7, 2.1655 /xm) 


5.25 ±0.55 


10.1 ±0.2 


9.04 ±0. 


.25 


148 ±6 


H 2 1-0 S(0) (2.2227 /xm) 


0.67 ±0.14 


0.11 ±0.03 


0.25 ±0. 


07 


9.8 ± 1.6 


H 2 2 - 1 S(l) (2.2471 /xm) 


0.29 ±0.10 


0.28 ±0.03 


0.20 ±0. 


06 


6.3 ± 1.2 


H 2 1-0Q(1) (2.4059 /xm) 


0.92 ±0.27 


1.04 ±0.22 


1.31 ±0. 


32 


13.2 ±3.2 


Equivalent width of stellar absorption features d 


(A) 






Si I (1.59 /xm) 


3.6 


3.3 


3.7 




3.4 


12 CO(6,3) (1.62 /xm) 


5.6 


3.4 


4.6 




4.8 


12 CO(2,0) (2.29 /xm) 


15.2 


8.4 


12.2 




13.0 



a Position of the 2.25" x 2.25" apertures defining the individual regions, relative to the nucleus. 
b Estimated uncertainties are 15% and 10% for the H- and if -band flux densities, respectively. 

c Quoted uncertainties represent those of the continuum subtraction (see section 3.1). Ionic tran- 
sitions are given as lower level - upper level. H 2 transitions are labeled by the upper and lower 
vibrational quantum numbers followed by S(j) or Q(j) which refer to transitions for which j — f 
equals —2 or respectively, where j and j' are the lower and upper rotational quantum numbers. 

Uncertainties on the equivalent widths are ±0.4 A for Si I, ±0.3 A for 12 CO (6,3), and ±0.6 A 
for 12 CO (2,0). The Si I equivalent widths are corrected for dilution by Brl4 (see section 3.1). 
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Table 3 

Comparison between various BR7 flux measurements 



Aperture diameter a 


Observed flux 


Reference 


(arcsec) 


(1CT 17 Wm- 2 ) 




3.8 


11.3 ±1.8 


1 




20.0 ±0.4 


2 




5.7 ±0.7 


3 




17.1 ±3.4 


4 


8 


48 ±8 


l b 




22 ±2 


5 




26 ±3 


3 




73.6 ± 14.7 


4 



a Aperture centered on the nucleus of M82. 

b An 8" aperture slightly exceeds the regions cov- 
ered by the 3D Br7 map, but we estimate the missing 
flux to be less than a few percents. 

References. — (1) This work; (2) Lester et al. 
1990; (3) Larkin et al. 1994; (4) Satyapal et al. 1995; 
(5) Rieke et al. 1980. 
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Table 4 

isosws line measurements in m 82 



Species Transition a A bsorvod FWHM Flux b Observation c Scaled flux d 







(Atm) 


(/mi) 


(10~ 15 Win" 2 ) 










(io- 15 W 


H 2 


1 - Q(3) (2.4238 fim) 


2.42539 


0.00214 


0.27 


02 


14" 


X 


20" 




H I 


Br/3 (n = 4- 6, 2.6252 M m) 


2.62726 


0.00295 


3.90 


01 


14" 


X 


20" 








2.62720 


0.00231 


4.10 


02 


14" 


X 


20" 




HI 


PM (n = 5 - 9, 3.2961 jim) 


3.29913 


0.00325 


0.59 


01 


14" 


X 


20" 




H I 


Pf 7 (n = 5 - 8, 3.7395 jum) 


3.74247 


0.00501 


1.07 


01 


14" 


X 


20" 




HI 


Hul4 (n = 6 - 14, 4.0198 /xm) 


4.02301 


0.00370 


0.13 


02 


14" 


X 


20" 




HI 


Bra (n = 4-5, 4.0512 /xm) 


4.05437 


0.00464 


8.75 


01 


14" 


X 


20" 








4.05474 


0.00414 


8.15 


02 


14" 


X 


20" 




H I 


Pf/3 (n = 5 - 7, 4.6525 Aim) 


4.65623 


0.00447 


1.39 


01 


14" 


X 


20" 




H 2 


0-0 5(7) (5.5112 jim) 


5.51437 


0.00657 


0.48 


02 


14" 


X 


20" 




H 2 


0-0 5(5) (6.9095 Mm) 


6.91299 


0.01155 


1.08 


01 


14" 


X 


20" 








6.91339 


0.00684 


1.15 


02 


14" 


X 


20" 




[Ar II] 


2 P 3 /2- 2 Pi/2 (6.9853 A*m) 


6.98912 


0.01083 


26.7 


01 


14" 


X 


20" 




H I 


Pfa (n = 5 -6, 7.4578 /an) 


7.46376 


0.00872 


2.49 


01 


14" 


X 


20" 








7.46417 


0.00727 


2.59 


02 


14" 


X 


20" 




H 2 


0- 5(4) (8.0251 Jim) 


8.02952 


0.00858 


0.68 


01 


14" 


X 


20" 




[Ar III] 


3 P 2 - 3 Pi (8.9914 Aim) 


8.99664 


0.01096 


4.89 


01 


14" 


X 


20" 








8.99681 


0.00824 


4.76 


02 


14" 


X 


20" 




H 2 


0- 5(3) (9.6649 Aim) 
2 Pi/2- 2 P 3 /2 (10.5105 Aim) 


9.67182 


0.01009 


0.97 


01 


14" 


X 


20" 




[S IV] 


10.51646 


0.01468 


1.89 


01 


14" 


X 


20" 






10.51639 


0.00825 


1.49 


02 


14" 


X 


20" 




H 2 


0- 5(2) (12.2786 Aim) 


12.29064 


0.02298 


2.00 


01 


14" 


X 


20" 








12.27736 


0.01598 


1.14 


02 


14" 


X 


20" 




[Ne II] 


2 P 3 /2- 2 Pi/2 (12.8136 Aim) 


12.81923 


0.01554 


99.1 


01 


14" 


X 


27" 


79.3 




12.81936 


0.01216 


89.2 


02 


14" 


X 


27" 


71.4 


[Ne III] 


3 P2 - 3 Pi (15.5551 Aim) 


15.56331 


0.01685 


17.5 


01 


14" 


X 


27" 


14.0 






15.56337 


0.01404 


15.7 


02 


14" 


X 


27" 


12.6 


H 2 


0- 5(1) (17.0348 Aim) 


17.04484 


0.01486 


1.38 


01 


14" 


X 


27" 








17.04508 


0.01228 


1.50 


02 


14" 


X 


27" 




[S III] 


3 Pi - 3 Pj (18.7130 Aim) 


18.72210 


0.01734 


34.7 


01 


14" 


X 


27" 


27.8 






18.72181 


0.01384 


31.5 


02 


14" 


X 


27" 


25.2 


[Ar III] 


3 Pi - 3 P (21.8293 Aim) 


21.84123 


0.05191 


0.70 


02 


14" 


X 


27" 


0.6 


H 2 


0- 5(0) (28.2188 Aim) 


28.23188 


0.03800 


0.78 


02 


20" 


X 


27" 




[S III] 


3 P - 3 Pi (33.4810 Aim) 


33.49818 


0.03873 


83.3 


01 


20" 


X 


33" 


58.3 






33.49816 


0.03432 


80.3 


02 


20" 


X 


33" 


56.2 


[Ne III] 


3 Pi - 3 P (36.0135 Aim) 


36.03283 


0.04118 


2.92 


01 


20" 


X 


33" 


2.04 






36.03075 


0.03535 


2.68 


02 


20" 


X 


33" 


1.88 



a Ionic transitions are given as lower level - upper level. H 2 transitions are labeled by the upper and lower vibrational 
quantum numbers followed by S(j), Q{j), or O(j) which refer, respectively, to transitions for which j — f equals —2, 
0, or 2, where j and j' are the lower and upper rotational quantum numbers. 

b The uncertainties on the observed line fluxes are estimated to be 20% for most of the lines, and up to 50% for 
the faintest lines (including uncertainties of the absolute calibration, continuum subtraction, and systematic errors). 

c 01: data from the full scan SWS01 spectrum, 02: data from the individual line scan SWS02 spectra. The aperture 
size is given as well. The SWS02 data are adopted foMthe analysis because of their higher S/N ratio. 

d H recombination and ionic fine-structure lines obtained in the 14" x 27" and 20" x 33" apertures are scaled down 
to the 14" x 20" aperture using the factors 0.8 and 0.7 respectively (see section 3.3). The uncertainties on the scaling 
factors are estimated to be 10%. 



Table 5 

Hydrogen recombination line fluxes in M 82 



Line 


Observed flux a 
(10~ 16 Wm- 2 ) 


Aperture b 


Scaled flux c 
(10~ 16 Wm" 2 ) 


Reference 


Ha 0.6563 Atm 


450 


90" x 90" 


170 ± 68 


1 


Ha 0.6563 


510 


64" 


200 ± 82 


2 


Pa/3 1.2818 /im 


86 ± 17 


30" 


43 ± 10 


3 


Brl3 1.6109 Aim 


1.32 ±0.50 


3D 


1.5 ±0.6 


4 


Rr1 2 1 6407 um 

I J L X _ 1 .Ut:U 1 LXj 111 


i oq + 70 


3D 


1.5 ± 0.8 


4 


Brll 1.6807 Aim 


1.43 ±0.50 


3D 


1.6 ±0.6 


4 


BrlO 1.7362 Aim 


3.22 ±0.70 


3D 


3.7 ±0.9 


4 


Br7 2.1655 Aim 


4.8 ±0.8 


8" 


17 ±4 


4 


Br/3 2.6252 Aim 


41.0 ± 10.5 


SWS 


41.0 ± 10.5 


4 


Pf<5 3.2961 Aim 


5.9 ±2.0 


SWS 


5.9 ±2.0 


4 


Pf 7 3.7395 Aim 


10.7 ±3.1 


SWS 


10.7 ±3.1 


4 


Hul4 4.0198 Aim 


1.3 ±0.4 


SWS 


1.3 ±0.4 


4 


Bra 4.0512 Aim 


81.5 ±14.7 


SWS 


81.5 ± 14.7 


4 


Pf/3 4.6525 Aim 


13.9 ±4.1 


SWS 


13.9 ±4.1 


4 


Pfa 7.4578 Aim 


25.9 ±4.7 


SWS 


25.9 ±4.7 


4 


H27a 0.95 mm 


(2.11 ±0.43) x 10~ 2 


Total 


(8.1 ±2.2) x 10~ 3 


5 


H30a 1.29 mm 


(3.20 ±0.38) x 10~ 3 


21" 


(2.0 ±0.4) x 10~ 3 


6 


H40a 3.03 mm 


(9.22±i"||) x IO" 4 


19" 


(6.2 ±1.6) x 10~ 4 


7 


H41a 3.26 mm 


(1.32 ±0.22) x 10" 3 


Total 


(5.1 ± 1.3) x 10~ 4 


5 


H53a 6.98 mm 


(1.87 ±0.28) x 10~ 4 


41" 


(8.6 ± 2.0) x 10~ 5 


7 



a Uncertainties include those of the absolute flux calibration, continuum subtraction, 
and systematic errors whenever possible. For the Ha fluxes, none were given in the 
references, but those of the beam size corrections should dominate for the scaled fluxes. 

b All apertures except the 3D and SWS fields of view are centered on the nucleus 
of M82. "Total" refers to measurements integrated over the entire emission regions in 
M82. 

c Line fluxes after beam-size correction to match the SWS 14" x 20" aperture (see 
section 4.1.1). The uncertainties account for those of the measurements and of the beam 
size correction. 

References. — (1) McCarthy et al. 1987; (2) Young et al. 1988; (3) Average of the 
results derived from the Pa/3 fluxes of McLeod et al. 1993 and Satyapal et al. 1995 (see 
section 4.1.1); (4) this work; (5) Seaquist et al. 1996; (6) Seaquist et al. 1994; (7) Puxley 
et al. 1989. 
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Table 6 

Hydrogen recombination line ratios and best-fit global extinction in M 82 



Ratio or property Intrinsic UFS, Draine MIX, Draine UFS, GC MIX, GC a 

Ay = 4 mag Ay = 43 mag Ay = 4 mag Ay = 52 mag 



Ha/Bra 


2.92 x 10 1 


3.59 x 10 1 


3.60 x 10 1 


3.39 x 10 1 


3.02 x IO 1 


Ha/Bra 


2.92 x 10 1 


4.23 x 10 1 


4.24 x 10 1 


3.99 x 10 1 


3.56 x IO 1 


Pa/3/Bra 


1.77 x 10° 


1.23 x 10° 


2.98 x 10° 


1.16 x 10° 


2.50 x 10° 


Brl3/Bra 


4.48 x 


10- 2 


3.10 x IO" 2 


6.97 x IO" 2 


2.93 x IO" 2 


5.84 x IO" 2 


Brl2/Bra 


5.72 x 


10- 2 


3.04 x IO" 2 


6.73 x IO" 2 


2.86 x IO" 2 


5.64 x IO" 2 


Brll/Bra 


7.48 x 


10-2 


3.16 x IO" 2 


6.89 x IO" 2 


2.98 x IO" 2 


5.78 x IO" 2 


BrlO/Bra 


1.01 x 


HT 1 


7.06 x IO" 2 


1.51 x IO" 1 


6.66 x IO" 2 


1.26 x IO" 1 


Br7/Bra 


3.17 x 


HT 1 


2.70 x IO" 1 


4.75 x IO" 1 


2.55 x IO" 1 


3.95 x IO" 1 


Br/3/Bra 


5.35 x 


HT 1 


5.83 x 10- 1 


8.51 x IO" 1 


5.51 x IO" 1 


6.99 x IO" 1 


Pft/Bra 


8.75 x 


10-2 


7.66 x IO" 2 


9.03 x IO" 2 


7.77 x IO" 2 


9.43 x IO" 2 


Pf7/Bra 


1.29 x 


io- 1 


1.34 x IO" 1 


1.42 x IO" 1 


1.34 x IO" 1 


1.41 x IO" 1 


Hul4/Bra 


1.45 x 


10-2 


1.60 x IO" 2 


1.61 x IO" 2 


1.60 x IO" 2 


1.61 x IO" 2 


Pf/3/Bra 


2.01 x 


io- 1 


1.66 x IO" 1 


1.51 x IO" 1 


1.71 x IO" 1 


1.72 x IO" 1 


Pfa/Bra 


3.33 x 


io- 1 


2.92 x IO" 1 


2.16 x IO" 1 


3.11 x IO" 1 


2.89 x IO" 1 


H27a/Bra 


2.56 x 


10- 5 


8.73 x 10- 5 


5.36 x IO" 5 


8.24 x IO" 5 


3.72 x 10- 5 


H30a/Bra 


1.33 x 


10- 5 


2.16 x 10- 5 


1.32 x IO" 5 


2.03 x IO" 5 


9.19 x IO" 6 


H40a/Bra 


2.48 x 


10~ 6 


6.68 x 10- 6 


4.10 x IO" 6 


6.31 x IO" 6 


2.85 x IO" 6 


H41a/Bra 


2.04 x 


io- 6 


5.50 x IO" 6 


3.38 x 10- 6 


5.19 x IO" 6 


2.34 x IO" 6 


H53a/Bra 


4.37 x 


io- 7 


9.27 x IO" 7 


5.69 x IO" 7 


8.75 x IO" 7 


3.95 x IO" 7 


F£ ra (Wm-*) 






9.28 x 10- ib 


1.51 x 10" i4 


9.83 x 10" ib 


2.18 x IO" 14 


lOg < Qly C (S- 1 ) > b 








53.701°;}° 


53.50^0.40 


ro 7 O+0.08 



Note. — The results reported are the line fluxes relative to Bra corrected for the best-fit ex- 
tinction. The intrinsic line fluxes can be recovered using the dereddened Bra flux (F# ra ). "UFS" 
and "MIX" stand for uniform foreground screen and mixed model, respectively. "Draine" indicates 
that the extinction law of Draine 1989 was assumed throughout the infrared regime, while "GC" 
indicates that is was amended by the Galactic Center law (Lutz 1999) between 3 fim and 10 /um. 

a Adopted extinction model and 3 — 10 /itm extinction law (section 4.1.3). 

b < Qhyc > i s tne average Lyman continuum photon emission rate derived from the extinction- 
corrected line fluxes. The uncertainties represent the dispersion of the individual values. 
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Table 7 

Extinction towards the ionized gas for selected regions in M 82 



Property 


Central 35 pc 


Bl 


B2 


3D field 


AUfs (mag) 
A™ IX (mag) 


10 ±5 
23 ± 10 


8±2 
27 ±7 


11 ±2 
45 ±20 


9±3 
36 ± 16 


log(0 Lyc [s- 1 ]) a 


52.03^;^ 


52.37±^ 


52.511q'23 


53.64±££ 



a Ql yc derived from the observed Br7 fluxes corrected for the adopted 
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Table 8 

Gas-phase abundances of Ne, Ar, and S in M 82 



Line 






h 




Abundance b 




(Wm- 2 ) 


(erg cm 3 s" 








Bra 


2.38 x 10- 


14 


2.29 x 10" 


-26 






[Ne II] 12.8 /zm 


1.37 x 10" 


13 


1.05 x 10" 


-21 


Ne+ /H = 


: 1.26 x 10~ 4 


[Nelll] 15.6 /xm 


2.17 x 10- 


14 


1.89 x 10" 


-21 


Nc++/H 


= 1.11 x 10" 5 


[Nelll] 36.0 


2.41 x 10- 


15 


4.71 x 10" 


-22 


Ne++ /H 


= 4.92 x 10~ 6 












[Ne/H] ft 


0.08 dex 


[Ar II] 6.99 yum 


6.51 x 10- 


14 


1.55 x 10" 


-20 


Ar+/H = 


: 4.03 x 10~ 6 


[Ar III] 8.99 /xm 


1.70 x 10- 


14 


1.25 x 10" 


-20 


Ar++/H 


= 1.31 x 10~ 6 


[Arlll] 21.8 yum 


1.09 x 10- 


15 


4.39 x 10" 


-21 


Ar++/H 


= 2.38 x 10" 7 












[Ar/H] « 


0.1 dex 


[S III] 18.7 Mm 


5.21 x 10- 


14 


1.87 x 10" 


-20 


S++/H = 


: 2.69 x 10~ b 


[S III] 33.5 yum 


7.34 x 10- 


14 


1.48 x 10" 


-20 


S++/H = 


= 4.79 x 10~ 6 


[S IV] 10.5 /xm 


5.70 x 10- 


15 


5.14 x 10" 


-20 


S+++/H 


= 1.07 x 10~ 7 












[S/H] ft - 


-0.6 dex 



a Fluxes for the 14" x 20" SWS aperture corrected for A^ lx = 52 mag. 
The Bra flux is computed from the Lyman continuum photon emission rate 
(table 6). 

b Ionic or elemental number abundance relative to hydrogen, with es- 
timated uncertainties of approximately ±50%, and up to a factor of two 
for the ionic abundances derived from the weakest lines. The total ele- 
mental abundance is given as [X/H] = log(X/H)M82 — log(X/H) Q , with 
Ne/H = 1.17 x 10~ 4 , Ar/H = 3.98 x 10~ 6 , and S/H = 1.62 x 10~ 5 for 
the solar composition (Grevesse & Anders 1989 and Grevesse & Noels 1993, 
as quoted by Ferland 1996). 



38 



Table 9 

Effective temperature and number of OB stars for selected regions in M 82 



Region 


Diagnostic 


Ratio a 


(K) 


iV 08 .5V b 


Starburst core c 






37400 


2.4 x 10 5 


SWS 14" x 20" 


[Nelll] 15.6/tm/[Ne II] 12.8 /im 


0.16 ±0.04 


37400 ± 400 


1.2 x 10 5 




[At III] 8.99 /im/ [Aril] 6.99 /an 


0.26 ±0.08 


33500 ± 500 






[S IV] 10.5//m/[S III] 18.7 /im 


0.11 ±0.04 


39900 ± 1000 




Central 35 pc 


He I 2.058/Br 7 


0.52 ±0.08 


35700 ± 800 


2040 




He I 1.701/BrlO 


< 0.13 


< 34200 




Bl 


He I 2.058/Br 7 


0.55 ±0.02 


36000 ± 200 


4470 




He I 1.701/BrlO 


0.22 ±0.06 


36400 ± 1600 




B2 


He I 2.058/Br 7 


0.52 ±0.02 


35600 ± 200 


6170 




He I 1.701/BrlO 


0.22 ±0.06 


36400 ± 1600 




3D field 


He I 2.058/Br 7 


0.55 ±0.03 


36000 ± 300 


8.3 x 10 4 




He I 1.701/BrlO 


0.20 ±0.09 


36000 ± 2300 





a The ratios are corrected for extinction and, when appropriate, for aperture size. The 
uncertainties include those of the continuum subtraction, systematic effects, and extinction 
correction. 

b A'os.sv is the number of equivalent 08.5 V stars required to produce the derived intrinsic 
Lyman continuum photon emission rates (see section 4.4.3). 

c The T g? and representative type for the SWS field of view is adopted for the entire 
starburst core (see section 4.4.3). 
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Table 10 

Properties derived from the near-infrared continuum emission for selected regions 



Property 


Units 


Central 35 pc 


Bl 


B2 


3D field 


Starburst core 


Spectral type 




K5I 


K2 I 


K4I 


K4I 


K4I 


£>1.6 


% 

















£>2.3 


% 





25 


15 


10 





D°h B 


% 





2 


1 








DT 


% 





1 


1 








D Ncb 


% 


< 1 


< 8 


< 5 


< 3 





riNcb 
K 


% 


< 1 


< 7 


< 4 


< 3 







mag 


10 ±4 


6±3 


8±3 


6±3 


9 


^MIX 


mag 


45 ±20 


25 ± 10 


28 ± 10 


17±7 


21 


L K 


10 8 L 


0.56 ±0.23 


0.074 ±0.024 


0.15 ±0.05 


2.7 ±0.8 


13 ±4 






1.1 x 10 4 


4400 


5800 


1.0 x 10 5 


5.0 x 10 5 



Note. — Di, 6 and D2.3 are the total amounts of dilution near 1.6 [im and 2.3 ^m, with typical 
uncertainties of ±10%. Z?g B and D K B arc the contributions to the H- and X-band flux densities 
from OB stars while D™ ch and D^ ob are those from nebular free-free and free-bound processes. 
The extinction values are those derived from the spectral fits except for the starburst core, which 
correspond to the effective extinction inferred as described in appendix B. We adopt the results for 
the mixed model throughout this work. Lk is the intrinsic K-hand luminosity from the evolved 
stars corrected for dilution and extinction while N* is the number of representative stars required to 
produce Lk- 
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Table 11 

Properties relevant in the determination of the effective ionization parameter 



Property 


Units Bl 


B2 


Starburst core 




Geometrical properties 




Sphere radius 


pc 19.5 


19.5 


225 


Column/disk radius 


pc 19.5 


19.5 


200 


Column/disk length 


pc 230 


365 


200 




Gas properties 






Mci 


M© 220 


220 


545 


r c \ 


pc 0.4 


0.4 


0.6 


Mh 2 


M© 1.1 x 10 6 


2.3 x 10 6 


1.8 x 10 8 


M H+ 


M© 1.4 x 10 5 


2.0 x 10 5 


7.6 x 10 6 


N cl 


5000 


10450 


3.3 x 10 5 


n c \ 


10" 2 pc" 3 16-1.8 


34 - 2.4 


0.69-1.3 




pc 2.3-4.7 


1.8-4.3 


6.5-5.3 


Ki 


cm- 2 5.5 x 10 22 


5.5 x 10 22 


6.0 x 10 22 


n 


pc 1.0 


0.88 


0.99 




Ionizing clusters properties 




^Lyc 


s- 1 2.34 x 10 52 


3.24 x 10 52 


1.23 x 10 54 


iV* 


2660 


3690 


1.4 x 10 5 


n± 


10~ 2 pc- 3 8.6-0.97 


12-0.85 


0.29-0.56 


d+-+ 


pc 2.8-5.8 


2.5-6.1 


8.7-7.0 




Structural and nebular properties 




\ogU a 


-1.2 


-1.1 


-1.7 


logoff a 


-2.3 


-2.4 


-2.4 



Note. — Two limiting cases are considered for the volume of each region: a 
sphere or an edge-on disk for the starburst core, and spheres or columns through 
the edge-on disk for Bl and B2. The ranges for the space density and mean 
separation of the gas clouds and stellar clusters correspond to these limiting cases. 

a The log U is computed from the direct application of equation (5) with the 
total intrinsic Lyman continuum photon emission rates and assuming spherical 
volumes. The log U e g is derived from equation (A3) as appropriate for the ran- 
domized distribution of clouds and clusters within the starburst regions of M82. 
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Table 12 

bolometric luminosity estimates for selected regions in m 82 



Region 




T evolved stars 


T tot 

-^bol 




(10 8 L ) 


(10 8 L ) 


(10 8 L ) 


Central 35 pc 


5.1 ± 1.8 


11 ±7 


18 ±9 


Bl 


5.7 ±2.0 


1.5 ±0.9 


8.9 ±3.1 


B2 


7.1 ±2.5 


2.9 ± 1.7 


12 ±4 


3D field 


130 ±50 


55 ±31 


220 ± 90 


Starburst core 


300 ± 30 


270 ± 80 


660 ± 120 



Note. — Lir is the infrared luminosity between 5 /zm 
and 300 fim and L b ^ lved stars is the bolomctric luminos- 
ity from the evolved stellar population. includes 
Lm, Ll v ° lvcd stars , and an additional 30% of L m (±10%) 
to account for light escaping in directions perpendicular 
to the galactic plane of M82. The uncertainties account 
for those on the absolute calibration of the data used to 
derive the quantities, on the conversion factors between 
the 12.4 /im flux density and Lm,, and between Lk and 
^bo° lvcd stars i on tne extinction, on the hot dust contri- 
bution to Lk , and on the fraction of escaping light. 
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Table 13 

Dynamical, gaseous, and stellar mass estimates in M82 



Region 


Mdyn 

(10 6 M ) 


M H2 
(10 6 M ) 


M H + a 
(10 6 M ) 


M* 
(10 6 M ) 


Central 35 pc 
Starburst core 


80 ±20 
800 ± 200 


0.821J-" 
180 ±50 


0.06613:^8 

7 fi+ 27 

'•°-5.4 


791^2 

fi in+270 
Dlu -250 



a The uncertainties include those on the intrinsic Lyman contin- 
uum luminosity, and on the electron density and temperature. 
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